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ABSTRACT 

Context. It was recently demonstrated that contact binaries occur in globular clusters (GCs) only immediately below turn-off point 
and in the region of blue straggler stars (BSs). In addition, observations indicate that at least a significant fraction of BSs in these 
clusters was formed by the binary mass-transfer mechanism. 

Aims. The aim of our present investigation is to obtain and analyze a set of evolutionary models of cool, close detached binaries with 
a low metal abundance, which are characteristic of GC. 

Methods. We computed the evolution of 975 models of initially detached, cool close binaries with different initial parameters. The 
models include mass exchange between components as well as mass and angular momentum loss due to the magnetized winds for 
very low-metallicity binaries with Z = 0.001. The models are interpreted in the context of existing data on contact binary and blue 
straggler members of GCs. The models are based on our recently developed code for evolutionary modeling of cool close and contact 
binaries. 

Results. The model parameters agree well with the observed positions of the GC contact binaries in the Hertzsprung-Russell diagram, 
and we achieve detailed parameters of several individual contact and near-contact binaries. Contact binaries in the lower part of the 
cluster main sequence are absent because there are no binaries with initial orbital periods shorter than 1.5 d. Contact binaries end 
their evolution as mergers that appear in the BS region. Binary-formed BSs populate the whole observed BS region in a GC, but a 
gap is visible between low-mass mergers that are concentrated along the zero-age main sequence and binary BSs occupying the red 
part of the BS region. Very few binary mergers are expected to rotate rapidly and/or possess chemical peculiarities resulting from 
the exposure of the layers processed by CNO nuclear reactions. All other binary mergers are indistinguishable from the collisionally 
formed mergers. The results show that binary-formed BSs may constitute at least a substantial fraction of all BSs in a GC. 
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1. Introduction 

W UMa-type stars are solar-type eclipsing variables with nearly 
sinusoidal light curves in which both binary components share 
a common envelope that is located between the inner and outer 
critical surfaces (equipotential zero-velocity surfaces crossing 
the LI and L 2 Lagrangian poi nts), also called inner and outer 
Roche lobes (lMochnackilll981l) . The orbital periods of W UM- 
type stars are concentrated around 0.25-0.6 d. Because of the 
physical properties of W UMa-types stars, we also call them 
contact binaries (CB). 

Numerous W UMa-type stars a r e obse rved in globular and 
old open clusters (iRucinskil Il998l l2000l) . but they are com¬ 
pletely absent among T Tauri stars and in young clusters, even 
though they are easy to detect (lMathieu|ll994 iMelo et ^l200lt 
iHebb et aDl2010l) . None is known in the intermediate-age clus¬ 
ter Hyades, in spite of t he high degree of completeness achieved 
in the binary search by iGriffinI (Il985h . and on ly one certain W 
UMa-type member is identified in Praesepe (IRucinskil Il998h . 
These data indicate that it takes several Gyr - up to the age of 
a globular cluster (GC) - before a CB is formed from its pro¬ 
genitor. The age of several Gyr als o results from the kinematic 
analysis of field W UMa-type stars dGuinan & Bradstreeljll988l: 
iBilir et al.ll2n05l) . 

It is generally accepted that progenitors of CBs are cool de¬ 
tached binaries. In the course of evolution they lose angular mo¬ 
mentum (AM) by a magnetized wind blowing from one or both 


components. With synchronous rotation, stellar AM loss (AML) 
results in a tightening of the orbit until the primary hlls its in¬ 
ner Roche lobe and Roche-lobe overflow (RLOF) occurs. For 
a binary whose initial orbital period is close to 2 d and whose 
initial primary is similar to the Sun, the time to reach RLOF is 
approximately equal to the e volutionary life time of the primary 
on the main sequence (MS) dSteniehlbOl lah . For shorter initial 
periods and/or lower primary mass, the MS life time is longer. 
Nonetheless, as long as the initial orbital period is not very dif¬ 
ferent from 2 d, the solar type initial primary is already advanced 
in evolution when RLOF occurs. 

W UMa-type stars observed in GCs always lie close to the 
turn-off (TO) in the Hertzsprung-Russell diagram (HRD), with a 
signihcant fraction falling above TO, in the region of blue strag¬ 
gler stars (BSs), wh ich are stars lyin g on the high-mass extension 
of the cluster MS dRucinskill2000l) . Early observations of GCs 
were not sufficiently sensitive to detect CBs in the lower part of 
the cluster MS, but the recent deep surveys of a few GCs revealed 
an absence of CBs fainter than about one magnitude below TO 
dKaluznv & ThomDsonI l200^ iKaluznv et al.ll20ld. l2013^ . As 
iKaluznv et al.ld2013all stressed, not a single system was detected 
in over 12 000 stars in the M4 cluster with magnitudes 17.8 
< V <21, that is, between one and four magnitudes below TO, 
whose light curves were analyzed. 

The aim of our present investigation is to obtain and ana¬ 
lyze a set of evolutionary models of cool, close, detached bi¬ 
naries with a low metal abundance, which is characteristic of 
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GC. Depending on the specific initial conditions, a fraction 
of the investigated binaries changes into CBs at the later age. 
Evolution ary computations are based on a model developed by 
one o f us (IStepiehl l 99511200^ iGazeas & Stepiehl2008l:[Stepied 
l2009ll . Altogether, about 10^ models haye been computed with 
an initial primary mass between 0.7-0.9 Mq, an initial secondary 
mass between 0.4-0. 8 Mq, an initial orbital period between 1.5- 
2.9 d, and a metal abundance Z = 0.001. The results are com¬ 
pared with the available observations of GCs, in particular with 
the data on the cluster CBs and BSs. 

The next section presents the essentials of the evolutionary 
model, including the recent improvements of the computer pro¬ 
gram, which follow the evolution of both components more ac¬ 
curately than in the original code (in particular, the slow mass- 
transfer rate is now calculated every time step), and additional 
effects are included. The third section describes the models to¬ 
gether with the statistical properties of the cool close binary pop¬ 
ulation as it evolves from the beginning to an age of 13.5 Gyr. 
Bulk properties of binaries that form mergers or are in a contact 
state at an age of 11, 12, and 13 Gyr are discussed. The ob¬ 
served properties of CBs and BSs, including individual objects 
from different GCs, are compared with the model predictions in 
the fourth section. The last section summarizes and discusses the 
results of the paper. 

2. Description of an evoiutionary model of cool 
close binaries 

2.1. Equations 

The cool close binary model (CCBM) describes the evolution of 
a cool close binary from the ZAMS untill a stage immediately 
preceding the merger of the components or the formation of an 
Algol-type star with an orbital period of several days. We do 
not follow the detailed evolution of the Algol-type binary here 
because it does not fulhll some of the CCBM assumptions. In 
addition, the Algol and post-Algol evolution has be en investi¬ 
gated thoroughly in the past and is well kno wn (e. g.. lPaczyhskil 
1 197 it ISarna|[l9^ iNelson & Eggleton]l200lh . 

As an initial close binary we took a detached system with 
an orbital period short enough to ensure full synchronization of 
the orbital and rotational periods. Because we are mostly inter¬ 
ested in progenitors of CBs, we discuss only binaries with initial 
periods shorter than 3 d. Longer period models always evolve 
into Algols. Cool means here that both components possess sub- 
photospheric convection zones that generate magnetic activity, 
in particular, a magnetized stellar wind. An approximate upper 
mass limit for these stars is about 1.3 Mq, although somewhat 
more massive stars can also be considered because they later de¬ 
velop a convection zone during their MS life. We here focus on 
stars with initial masses equal to or lower than 0.9 Mq. They 
have a sufficiently deep convection zone already on the ZAMS. 
We apply the Roche model to describe the orbital parameters. 

The basic equations of CCBM are the third Kepler law, the 
expression for AM, and t he approximate e xpressions for inner 
Roche-lobe sizes rj and ra lEggletorJ (Il983h . 



(1) 

^tot “ ^spin + ^orb ■' 

(2) 

where 


//,pi„ = 7.13 X lQ^°{k\MiR\ + klM2Rl)P ^^, 

(3) 


and 

//orb = 1-24 X + q)-^ , (4) 

a Q.6q^i^ + ln(l + q^!^) ’ 

^ ^ O.A9q-^l^ 

a 0.6q^^i^ + ln(l -H 

Here P is the period (rotational and orbital), M - Mi + M 2 is 
the total mass with component masses Mi and M 2 , /?iand R2 are 
component radii, a is the semi-axis, //tot, //spin and // 01 b are total, 
rotational and orbital AM, kj and are the (nondimensional) 
gyration radii of the two components, ri and r2 are the sizes of 
the inner Roche lobes, and q - Mi I M2 is the mass ratio. Masses, 
radii, and semi-axis are given in solar units, the period is given 
in days, and the AM in cgs units. 

We assumed that stellar winds from the two components are 
the dominating mechanism of the orbit evolution. They carry 
away mass and AM according to the formulas 

Mi,2 = -10-"/?2 2, (7) 

ATT 

= -4.9 X 10''\r\Mi + rIM2)IP . (8) 

Here M is in solar mass per year and d//tot/df is in gcm^s_' 
per yea r. The formulas were d i scusse d at length by ISteni^ 
(l2006allbl) and iGazeas & Steniehl (l2008h . They are calibrated by 
the observational data of the rotation of single, magnetically ac¬ 
tive stars and empirically determined mass-loss rates of single, 
solar type stars. Both formulas apply in a limiting case of a 
rapidly rotating star in the saturated regime. Note that they do not 
contain any free adjustable parameters. The constant in Eq. O 
is uncertai n within a factor of 2 and that in E q. ([8]l is uncertain 
to + 30% (ISteDiehll2006bt IWood et ani2002h . Any interaction 
between winds from the two components is neglected. 

2.2. Initial assumptions 

The model equations were applied to binaries with parameters 
characteristic of GC members. The CCBM evolution was fol¬ 
lowed until the age of the oldest GCs and snapshots of the model 
population properties are discussed in detail for three speci- 
hed ages of 11, 12 and 13 Gyr. We focused on stars with ini¬ 
tial masses close to the TO mass, because more massive clus¬ 
ter members have already completed their evolution and have 
formed compact objects, whereas the less massive ones are still 
on the MS, burning hydrogen in their cores. Specihcally, we took 
initial primary masses from the interval 0.7-0.9 Mq, with the de¬ 
tailed values given in Table [1] If the primary in a model has a 
mass that is closely spaced between 0.86-0.87 Mq , it means 
that the primaries reach RLOE between 11 and 13 Gyr. The ini¬ 
tial masses of the secondary component were adjusted to avoid 
extreme mass ratios qijin - 1 and qiait » 1. The evolution of bi¬ 
naries with exactly equal mass components is uninteresting to us, 
whereas the rapid mass transfer following RLOE is very likely 
nonconservative for high values of the mass ratio. Without a reli¬ 
able theory describing this process, the stand ard procedure is to 
use fr ee parameters for mass and AM losses (ISarna & De Grev3 
Il996h . To avoid introducing them, we assumed conservative 
mass transfer. Our models have 1.08 < qi^n < 2.25. The masses 
of the secondary component are also given in Table [1] Eor each 
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Table 1. Initial component masses of the computed binary mod¬ 
els. For each combination of Mi and M 2 a set of 15 evolution¬ 
ary models has been calculated with initial orbital periods 1.5, 
1.6,...,2.8,2.9d 


Mi(Mo) 

MiiMo) 

0.7 

0 . 4 , 0 . 5 , 0.6 

0.75 

0 . 4 , 0 . 5 , 0.6 

0.8 

0 . 4 , 0 . 5 , 0 . 6 , 0.7 

0.81 

0 . 4 , 0 . 5 , 0 . 6 , 0.7 

0.82 

0 . 4 , 0 . 5 , 0 . 6 , 0.7 

0.83 

0 . 4 , 0 . 5 , 0 . 6 , 0.7 

0.84 

0 . 4 , 0 . 5 , 0 . 6 , 0.7 

0.85 

0 . 4 , 0 . 5 , 0 . 6 , 0.7 

0.86 

0 . 4 , 0 . 5 , 0 . 6 , 0.7 

0.865 

0 . 4 , 0 . 5 , 0 . 6 , 0 . 7 , 0.8 

0.8675 

0 . 4 , 0 . 5 , 0 . 6 , 0 . 7 , 0.8 

0.87 

0 . 4 , 0 . 5 , 0 . 6 , 0 . 7 , 0.8 

0.88 

0 . 4 , 0 . 5 , 0 . 6 , 0 . 7 , 0.8 

0.89 

0 . 4 , 0 . 5 , 0 . 6 , 0 . 7 , 0.8 

0.9 

0 . 4 , 0 . 5 , 0 . 6 , 0 . 7 , 0.8 


pair of the component masses, 15 models were computed with 
initial period lengths from 1.5 d to 2.9 d, every 0.1 d. Altogether, 
975 evolutionary models were obtained. In the following, we 
identify individual models by giving the initial values of three 
basic parameters: Mi -i- M 2 (P), for instance, 0.8-i-0.6(1.8). 

We adopted the range of initial period values according to 
the following considerations. It is now generally accepted that 


binary and multiple stars are for med during the_earK 


frag - 

1994 


tne early 

mentation of a protostellar cloud dBossI 1 19931 : iBonnelll [l. 
iKratter et al.l 1201(1 iMachida et al.l l2008h . After the end of the 
rapid accretion phase, cool protostars enter the T Tauri phase 
when the mass has almost reached its final value and changes 
only li ttle. This transition takes place at the age of 10^ - 10^ 
years (IMachida et al.ll20iol: IVorobvo\ll2010l : lOhtani & Tsuribel 
I 2 OI 3 I) . The orbit of a binary with two freshly formed T Tauri- 
type components must be wide enough to accommodate both 
stars. The youngest T Tauri stars with masses bet ween 0.5 and 
1.5 Mp have radii from 2 up to over 3 solar radii dBaraffe et al.l 
ll998l:lTbgnelli et al.ll201ll) . which means that the initial orbital 
period must be longer than about 1.5 d (the accurate value de¬ 
pends only weakly on the mass ratio). We adopted this value as 
a lower limit for the initial orbital period of an isolated binary. 
Shorter period binaries can be formed by dynamical interactions 
(col lisions ) with other objects in a GC dHeggidl975l:lHillsll97'^ 
[Huil 19831 Ikeigh et alJlTOTl iHvoki & Giers7ll2ni 3h . but we ex¬ 
cluded them from our analysis. 

There exists another efficient mechanism that effectively 
shortens the orbital periods of some binaries. It requires a third 
body on a very wide and strongly inclined orbit around the in¬ 
ner binary. The third companion induces so-called Kozai cy¬ 
cles that periodically change the eccentricity of the inner or¬ 
bit. Tidal interactions between the components during perias- 
tron passages dissipate the system energy and tighten the orbit. 
Numerical models of this mechanism, called Kozai cycles with 
tidal friction (KCTF), show that a signihcant fraction of bina¬ 
ries with initial periods of 2-3 w eeks shorten their periods to 2- 
3 d w i thin several million ye^s (j ^gleton & Kiseleva-Eggletod 
20061: Fabrvckv & Tremain^ I^OTt Perets & Fabrvckvn2009t 
Naoz & Fabrvckvl 2014 . At thes e periods the orbits circular - 


ize and do not change any more dFabrvckv & Tremainell200% . 
Population synthesis indicates that among binaries with periods 
shorter than 3 d, those formed by the KCTF mechanism outnum¬ 


ber isolated systems (iFabrvckv & Trem ainel200' ^. Observations 
of field binaries conhrm this expectation. Tokovinin et al1(l2006ll 
found out that 96% of binaries with periods shorter than 3 d 
possess a third companion, in contrast to 34% of binaries with 
periods longer than 1 2 d. A similar conclusion was reached by 
iRucinski et al.l (l2007h . who showed that all or almost all held W 
UMa-type binaries possess a tertiary component. 

Consequently, we assumed that the initial period distribution 
of cool close binaries has a sharp maximum at 2-3 d, resulting 
mostly from the KCTF mechanism, and a rapid decline toward 
1.5 d. The extension of the distribution beyond 3 d is unimpor¬ 
tant for our purposes because, as we shown below, these binaries 
do not form CBs. 


2.3. Description of evoiutionary computations 

The ev olution of cool close binaries can be divided into three 
phase s (IStepiehl2006al:l^czvhski et al.ll200'^ iGazeas & Stepie"^ 
l2008h . where phase I begins at the initial state to RLOF by the 
initial primary (denoted with the subscript “1” throughout, also 
when it becomes a less massive component). The other com¬ 
ponent is marked “2”. The binary is detached in phase I, and 
its evolution is dominated by two mechanisms: individual evo¬ 
lution of each component across the MS, and the magnetized 
winds carrying away mass and AM. Any interaction between 
winds from the two components is neglected, and no other mech¬ 
anism influencing the orbit is considered. Equations ([T]l8]l are in¬ 
tegrated at every time step. The typical length of the time step 
was set to 6.8 x 10^ years so that the age of a GC is reached 
in about 200 steps. At each time step, the radii of the compo¬ 
nents were calculated using the grid o f single-star models w ith a 
metal content Z = 0.001 published bv lGir^di et al.l (l2000l) . The 
value of the radius was interpolated in time and mass. The to¬ 
tal mass lost during phase I amounts to several percent of the 
initial stellar mass, depending on stellar radius and the duration 
of phase I. This has a minor effect on the radius evolution of 
the components. More important is the evolutionary advance¬ 
ment due to hydrogen burning, which slowly inflates a star. At 
the same time, the Roche lobes descend onto the stellar surfaces 
due to AML (Eq. |8]l. If RLOE did not occur within 13.5 Gyr, 
calculations were stopped and the binary remained detached. To 
allow fo r the supersaturation of magnetic activity in ultra-fast 
rotato rs dRandich et al.lll995 iProsser et al.lll996l : IStepieri et al.l 
I 2 OOII) . the period in Eq. ([^i was replaced by 0.4 d for binaries 
with orbital periods shorter than that limiting value. Phase I is 
the longest evolutionary phase with a duration of the same order 
as the MS lifetime of the initial primary. 

Eollowing RLOE, phase II begins when a rapid mass trans¬ 
fer occurs from the initial primary (henceforth loser) to the initial 
secondary (gainer). It is much shorter than other phases and lasts 
no more than one or a fe w percent of phase I. The d etails of this 
process are discussed bv IStepieh & Kiragal (l2013ah . The gainer 
swells after receiving a small amount of matter, until it hlls its 
inner Roche l obe and the syste m temporarily assumes a contact 
conhguration (IWebbinklll976h . The mass that then flows from 
the loser encircles the s wollen gainer c lose to its equator and re¬ 
turns to the parent star (IStepiehll2009l) . With no mass accreted, 
the swollen star slowly shrinks, trying to revert to thermal equi¬ 
librium. When its surface moves below the Roche lobe, the star 
can again receive a small amount of mass, which causes it to ex¬ 
pand back to the size of the lobe. The process of shrinking, gain¬ 
ing some mass, and expanding again is repeated (a sort of equi¬ 
librium is established where an instantaneous mass-accretion 
rate just balances swelling and shrinking of the star close to the 
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Fig. 1. HR diagram with representative initial binary models 
plotted together with a low-metal isochrone assumed to track 
theZAMS ofaGC. 


Roche lobe), resulting in a mass-transfer rate governed by the 
thermal timescale of the gainer rather than loser. As soon as the 
thermal timescale of the loser becomes longer than that of gainer 
(this occurs when the two components have approximately equal 
masses), the gainer can receive all the transferred mass. Phase II 
ends when both stars regain equilibrium and the loser just fills 
its Roche lobe. The mass ratio is inverted so that the initial pri¬ 
mary is now less massive than the initial secondary. Regardless 
of the details of the mass transfer process, we assumed that mass 
is transferred on a thermal timescale of ~ 10® years. A constant 
mass transfer rate Mq = 5 x IO^^Mq yr“' was adopted based on 
observation that the total mass transferred during phase II in any 
of the modeled binary did not exceed 0.5 Mq. 

In most cases, a short-period semi-detached configuration 
of the Algol-type was obtained at the end of phase II with the 
gainer inside its Roche lobe. We call this a near contact binary 
(NCB). But a few binaries were so compact that a CB was al¬ 
ready formed before phase III began. 

Similarly as in phase I, the evolutionary expansion of both 
stars, accompanied by mass and AM loss by the winds, drives the 
binary evolution in phase III. In particular, the expansion of the 
loser beyond its Roche lobe results in a slow mass transfer (on 
the evolutionary timescale) to the gainer. This widens the orbit, 
which in turn suppresses the transfer. In absence of winds, the 
binary period systematically increases until mass transfer stops 
completely. In the presence of a low AML rate, the period still in¬ 
creases, although at a lower pace. If the AML rate is high and/or 
the mass transfer rate low, however, the orbit shrinks until both 
components merge. Each binary follows one of these routes, de¬ 
pending on the relative efficiency of the two processes. 

To compute the binary evolution in phase III, we used the 
same equations and assumptions as in phase I. Particularly im¬ 
portant was a proper treatment of the loser. Because we lack de¬ 
tailed evolutionary models of binary components after a rapid 
mass exchange, single-star models were used, as described 


above. Analytic formulas were fitted to mass-radius relations of 
stars at different evolutionary stages: at the MS, at the TAMS, 
and for several values of the helium core mass. In addition, for¬ 
mulas for time derivatives of stellar radii were found for differ¬ 
ent masses and evolutionary stages. Using them, the radius of 
the secondary component was computed at each time step and 
compared to the Roche lobe size. The mass transferrate resulted 
from this comparison. 

Several different formulas for Mtr following RLOF in cool 
close binaries hav e been used in the pas t assuming conserva¬ 
tive mass transfer. Tlben & Tutukov ( 1984t) assumed proportion - 
ality of Mtr to (Rt /ril^MEggleton & Ki seleva-EggletonI (l2002h 
to (ln(Ri/ri))®, whereas lifen et al.ril2000t) used 


Mtf = for Ri>n, (9) 

Mti = 0 for Ri < ri, 

with C - 1000 MQyr^' and a - 3. After several trials, we 
adopted the formula of Han et al., but with C = 1 MQyr ' and a 
- 5. The adopted values secured stability for most of the mod¬ 
els. Only in a few cases unnecessary oscillations appeared (see 
below). 

This method of calculating the mass-transfer rate was ap¬ 
plied to semi-detached binaries until the gainer reached its 
Roche lobe, thus producing a CB. After that, another approach 
had to be used because both stars overfill their Roche lobes. 
Now, the degree of overfill by the gainer was calculated at ev¬ 
ery time step and compared with the same quantity for the loser. 
If the former was lower than the latter, the mass transfer rate 
was increased according to Eq. 9, except that an upper limit of 
2 X 10“® MQyr^' was imposed to keep the changes of the binary 
configuration smooth. If not, the zero transfer rate was applied 
for the current time step. The model calculations were stopped 
when the degree of the inner Roche-lobe overfill by both stars 
reached 20%, attaining thus the outer Roche lobe for typical val¬ 
ues of the mass ratio of the modeled binaries {q » 0.6 - 0.8). 
At this moment, the binary was assumed to lose mass and AM 
through the Lagrangian L 2 point at a rapidly increasing rate 
with an ultimate merging of the components. Similarly, merg¬ 
ing should also occur as a result of the Darwin instability when 
the spin AM exceeds o ne third of th e orbital AM. This occurs for 
mass ratios below 0.1 (lRasiolll99^ . We stopped the calculations 
at this value, notwithstanding the low degree of overfill. 

The contact configuration arises from the NCB formed at the 
beginning of phase III when the influence of AML on the orbit 
prevails, or at least balances the mass-transfer rate so that the 
orbit remains tight during phase III. In the opposite situation, 
the orbit widens, the period increases, and the binary remains 
in the Algol configuration with the gainer well inside its Roche 
lobe. This occurs when the initial primary reaches RLOF after 
completing its MS evolution, so that hydrogen in the central part 
of its core is already exhausted. Consequently, it moves to the 
red giant region on the HRD, and its rapid expansion results in a 
high enough mass-transfer rate to dominate the orbit evolution. 
We stopped the calculations for orbital periods exceeding 3 d 
when the spin-orbit synchronization may already break down. 

The spin AM of the two components is about two orders of 
magnitude lower than the orbital AM for all considered binaries 
during phases I-II and for almost all binaries during phase III. 
Nonetheless, they were calculated at each time step. Constant 
values of gyration radii were adopted: - 0.10, which corre¬ 

sponds to a star with core hydrogen d epleted, and - 0.06 , 
corresponding to an unevolved MS star dRutten & Pvlvserll 19881: 
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IClaret & Gitnenejll992h . The value of is overestimated for 
early evolutionary phases, but it better reflects the internal struc¬ 
ture of a loser in phase III. In addition, the total binary luminosity 
and the average effective temperature of CCBMs were calculated 
at each time step. The sum of the component luminosities was 
adopted as the binary luminosity, and the binary effective tem¬ 
perature was obtained from the surface-weighted mean bright¬ 
ness of both stars. 


Lb = Ti+T2, (10) 

Tts:,^UI^n(T(R\ + Rl), ( 11 ) 

where Lb is the binary luminosity and Leff_b is the binary ef¬ 
fective temperature. 

Figure[T]shows the region of the HRD that is occupied by the 
initial CCBMs. To avoid crowding, only 19 models were plot¬ 
ted with a primary mass between 0.7-0.9 Mq, every 0.1 Mq, 
and all secondary masses accompanying them. We also plot the 
isochrone f = 10^ years, Z = 0.001 from the latest Padova mod¬ 
els (http://stev.oapd.inaf.it/cmd), which shows the approximate 
position of the ZAMS for a GC. The modeled binaries all lie 
above the cluster MS, as expected. 


3. Results 

3.1. Period evolution 

The binary evolution sensitively depends sensitively on whether 
the primary mass is higher or lower than the TO mass. Primaries 
with lower masses remain detached untill the specified age of a 
GC, unless their initial orbital period is significantly shorter than 
1.5 d. Figure |2] presents an example of the period evolution of 
binaries with a primary mass lower than the TO mass. Out of all 
60 models with an initial primary mass of 0.8 Mq , 19 (32%) 
reach phase III within an age of 13.5 Gyr, whereas the others 
remain detached. Corresponding figures for other low masses 
look similar, except that the number of models reaching phase 
III increases from only 4 (9% ) for 0.7 Mq to 35 (58%) for 0.86 
Mq (note that because of the mass lost through stellar winds, 
the final masses of the stars are about 6-7% lower than initial). 
All binaries with initial primaries more massive than 0.86 Mq 
reach RLOF within the age of 13.5 Gyr and show pronounced 
evolutionary effects (see an example in Fig. [3]). Flere we also see 
that binaries with periods longer than ~ 2 d reach RLOF when 
their primaries have already left the MS, which results in a high 
mass-transfer rate in phase III and a rapid period increase. Some 
models with q approaching 0.1 show spurious period oscillations 
before merging, see for example the model in the lower right of 
Fig. [3] This is caused by a too long time step. We used a constant 
value of 5 X 10^ years for all models in phase III, which means 
that the typical duration of this phase is 200-400 steps. A high in¬ 
stantaneous mass-transfer rate, multiplied by the long time step, 
results in a rapid decrease of an already low loser mass, pro¬ 
ducing a sudden period jump, followed by its decrease at the 
next time step as the mass transfer rate drops to zero. We ob¬ 
tained exemplary calculations for the time step shortened four 
times, that is, for 1.25 x 10^ years. The period oscillations are 
then completely absent (Fig.Hl). Switching to such a short time 
step would, however, substantially lengthen the computations, 
therefore we decided to ignore these oscillations in the follow¬ 
ing analysis. We did so because a correct value of a binary pe¬ 
riod matters only when summation over consecutive bins is per¬ 
formed to obtain the period distribution at a specified age (see 
below). Instead of trying to determine a correct period value, for 



Fig. 2. Variation in time of the orbital period of all binary models 
with an initial primary mass equal to 0.8 Mq. Period variations 
of binaries that have not yet reached RLOF within the age of the 
Universe are shown with dotted lines and solid lines describe the 
period variations of binaries that reach RLOF within this age. 


instance, by determining the local minimum in the vicinity of 
this age, we included its instantaneous value in the summation. 
In effect, a small number of periods fall incorrectly into an adja¬ 
cent bin. This changes the general appearance of the distribution 
only little. The components merge as a result of the Darwin in¬ 
stability always near the local period minimum, which is close 
to the correct value (Fig.|4]i. 


3.2. Period distribution 

The orbital period distribution of cool close binaries in GC 
evolves in time under the influence of several mechanisms. 
Some, such as encounters with other stars, harden the bina- 


ries by shortening their period s (Heggie 1975; iHillsI 119751 : 
[Hu 3[T^ [Hvoki & Gierszfl^l iLeigh et al.ll2013h . or like the 


KCTF mechanism that operates on a timescale of several Gyr 
because of the very long i nitial period of the inner binary 
(iFabrvckv & Tremainel 1200^ . inject binaries into the consid¬ 
ered period interval. Other mechanisms such as expelling bina¬ 
ries from GC, component merging, or formation of Algols with 
longer periods, remove binaries from that interval. A full discus¬ 
sion of all these mechanisms is beyond the scope of the present 
paper. We consider here only the evolution of the period distri¬ 
bution of CCBMs. 

We made the following assumptions about the population of 
cool close binaries in a GC. 


1. The initial period distribution is constant for the period in¬ 
terval 2.0-2.9 d and decreases linearly to zero at 1.4 d: 

niP,Mi,q)dP = 0 for P < 1.4, 

n{PMi,q)dP = ^(P-L4) for L4<P<2, (12) 

n{PMi,q)dP = 1 for P > 2, 
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Fig. 3. Variation in time of the orbital period of all binary models 
with an initial primary mass equal to 0.88 Mq. All these binaries 
reached RLOF within the age of 11 Gyr. 



time(Gyr) 


Fig. 4. Period variations of a model 0.88+0.4(2.2) in phase III, 
calculated with a time step of 5 x 10® years (dotted line) and 
1.25 X 10® years, i. e. four times shorter (solid line). 


where n{P, M \, q) is the number density of binaries with dif¬ 
ferent periods and mass ratios. The value of 1.4 results from 
a linear approximation of the step-wise distribution consid¬ 
ered here, where a first non-zero value of n is for P = 1.5 
d. The lack of shorter initial periods is due to the inflated 
component radii during the T Tauri phase (see Sect. 2.2). 



P (day) 

orb 


Fig. 5. Evolution of the orbital period distribution of cool close 
binaries in time. Distributions are shown every 1 Gyr from 0 till 
13 Gyr, and they are supplemented with the distribution at 13.5 
Gyr. The ordinate gives the age of the binary population in Gyr. 


2. The initial period and component masses are not correlated, 
that is, n{P,M\,q)dPdMdq = ni(P)n 2 (Mi)n^(q)dPdMdq, 
where n\{P),n 2 {Mi) and n^iq) are the number densities de¬ 
pending on period, mass, and mass ratio, respectively. This 
follows from the semi-emp irically determ ined properties of 
the birth binary population (lKrouDalll995h . 

3. The primary mas s distribu t ion is given by the initial 
mass function of ISalneteil (Il955h : n 2 {M\) oc 
Al ternative forms of this func tion have als o been proposed 
by Miller & Scald (Il979t) and IScalol (Il986ll . More recently, 
iKrounal ^ 2002h confirmed the Salpeter mass function for stars 
with masses higher than 0.5 Mq , except that he suggested a 
somewhat flatter mass dependence with an exponent of 2.3. 
Nevertheless, the original Salpeter relation is still in use. The 
difference between both exponents has a negligible influence 
on the mass function over the narrow mass interval consid¬ 
ered here. 

4. The mass ratio distribution is flat, «r(q) = const. 
(iRaghavan et al.ll2010t iMilone et al.ll2012h . although short- 
period field binaries may show some preference to g' 1 
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dPuchene & Krausll2Q13h . In lack of detailed information on 

the mass ratio distribution of GC binaries, we adopted the 

flat distribution. 

Before calculating the period distribution, a weight for each 
CCBM must be determined to be used in the binary summation. 
The initial CCBMs are distributed uniformly in P, which means 
that each model represents binaries from an equal period inter¬ 
val. When summation is performed over consecutive period bins, 
the weight wp, attached to each model, is therefore constant for 
P > 2 d and decreases linearly to zero, as given by Eq. 12. The 
situation is more complicated when we consider weights wm and 
Wq connected with the initial primary mass and initial mass ra¬ 
tio, because the models are distributed nonuniformly in M\ and 
q. To obtain their values, it was necessary first to calculate the 
widths of the respective intervals in Mi and q, centered on each 
model, and then to multiply them by factors resulting from as¬ 
sumptions 3 and 4 above. Finally, a product of all three weights 
was attached to each of the 975 CCBMs and kept constant over 
all evolutionary calculations. 

Figure|5]presents the time evolution of the period distribution 
with the 0.1 d bin size. The evolving distribution is shown every 
1 Gyr, starting from the initial distribution at the bottom of the 
figure untill the age of 13 Gyr. The last distribution, at the top 
of the figure, corresponds to the age of 13.5 Gyr. Broken lines 
show zero lines of the consecutive distributions. 

The main evolutionary effects of the period distribution con¬ 
sist of the shift of the short-period limit toward the shorter peri¬ 
ods, together with a stretching or flattening of the distribution, 
except for its long-period end, where a deficit of binaries is vis¬ 
ible. The deficit is due to the period increase of several binaries 
beyond the three-day limit, which is not compensated for by the 
possible period shortening of binaries with initial periods > 3 d. 
They were not included in the present models, but we can esti¬ 
mate their approximate influence on the plotted distribution. All 
binaries with initial periods 3-3.5 d shorten them below 3 d dur¬ 
ing phase I, but those with primaries more massive than the TO 
mass quickly increase their periods above 3 d past RFOF (simi¬ 
larly as shown in Fig.|3]l. Fower mass binaries with the same ini¬ 
tial periods remain detached over the Hubble time and shorten 
their periods continuously, but lose only about 15-20% of the 
initial AM because the AMF rate is inversely proportional to the 
period, see Eq. |8] Their final periods fall between 2.4-3 d, so 
they would partly All in the deficit if they are included in the dis¬ 
tribution. Nontheless, a comparison between the observed and 
predicted distribution should be limited to periods shorter than 
2.4 d. Binaries with initial periods above 3.5 never enter the pe¬ 
riod interval plotted in Fig. |5] This means that the distribution 
does not go to zero at the long-period boundary, but its exact 
shape there depends on the initial period distribution beyond 3 
d, which is not considered. 

The general properties of the period distribution are insen¬ 
sitive to the bin size, for example, increasing the bin size twice 
results in a less ragged shape, but the broadening toward short 
periods and flattening of the distribution are still clearly visible. 

The short-period limit crosses the 1 d mark at the age of 4 
Gyr and 0.5 d mark at about 7 Gyr. This is a very interesting 
result. It shows that the value of the short-period limit of de¬ 
tached MS binaries in a stellar cluster can be used as an inde¬ 
pendent measure of the cluster age, similarly as t he rotation pe¬ 
riod of cool MS star s is a measure of stellar age (iBarne ^ l200^ 
iMeibom et al.ll20flT) . At present, very few clusters are known 
with a sufficiently complete sample of binary stars to safely de¬ 
termine a short-period limit. However, the increasing accuracy 


of photometric observations in the future will enable detecting 
light variations caused by the ellipticity of noneclipsing binary 
components. As a result, a high percentage of all close bina¬ 
ries will be detected in several clusters and a short-period limit 
of the period distribution can be accurately determined. Since 
close binaries can also be formed by collisions with other cluster 
members (ignored here), their possible effects on this limit will 
have to be quantified before the method is applied to a particular 
cluster. 

The flattening of the distribution with time is accompanied 
by a decrease of the total number of binaries within the consid¬ 
ered period interval that is due to two mechanisms: merging of 
the two components, and a period increase beyond 3 d of Algol- 
type binaries. As a result, the total weighted number of CCBMs 
(equivalent to the integral of the distribution over period) in the 
period interval 0.15-3.0 d decreases by about one third (more 
accurately, by 32.9%) at 13.5 Gyr, compared to the initial data. 
Unweighted number of models included in the distribution de¬ 
creased at that time from an initial 975 binaries to 398 (41%). 

We conclude that the predicted shape of the period distribu¬ 
tion of cool close binaries at the age of a GC is flat over the 
interval 0.1-2.5 d and bears no similarity to the initial distribu¬ 
tion. 

3.3. Mergers and contact binaries in a giobuiar duster 

As we showed above, the binary evolution past RLOF depends 
on the relative importance of slow mass transfer that length¬ 
ens the orbital period and of the AMF that shortens it. Binaries 
with rapidly increasing period remain in an Algol configura¬ 
tion, whereas the rest experience a contact phase with the ul¬ 
timate merging of the components if the timescale for a merger 
is shorter than the cluster age. 

The model calculations show that coalescence typically oc¬ 
curs in less than 1 Gyr after RFOF, and the typical lifetime of a 
merger in the region of BS is also of the order of 1 Gyr. We there¬ 
fore focused on binaries with a RFOF occurring at, or later than, 
10 Gyr, which then evolve toward a CB or NCB configuration. 
Binaries evolving to the Algol configuration are not discussed. 

We selected model binaries in which RFOF occurred be¬ 
tween 10 and 11 Gyr and which have orbital periods shorter than 
1 d at the age of 11 Gyr. In many of them (particularly those with 
the shortest periods), both components already All their Roche 
lobes, but in others only the loser does, whereas the gainer is 
within its lobe. The computed stellar radii result strictly from the 
evolutionary models without allowing for any effects that might 
influence them (e. g., magnetic fields or deviations from thermal 
equilibrium). The observed radii of cool, rapidly rotating stars 
are often larger than the model radii, so several real stars with 
the same global parameters as CCBMs may in fact be already 
CBs even if they look like NCBs in our simulations. We marked 
binaries with P < 0.5 d as CBs and those with 1 > P > 0.5 d as 
NCBs and indicate them with asterisks and plus signs, respec¬ 
tively. Numerical data on these models are given in Table |2] with 
the consecutive columns giving the name of the model, compo¬ 
nent masses, mass ratio, period, luminosity, effective tempera¬ 
ture, and remarks: CB - contact binary, NCB - near-contact bi¬ 
nary, MS - binary lies below TO, BS - binary lies in the BS 
region. The same data are plotted in Fig. |6] together with the 
ZAMS and all Gyr isochrone. In addition, the red boundary of 
the BS region we adopted is shown as the long-dashed line. This 
is the ZAMS line shifted upward by AlogL = 0.8 (= 2 mag). 
It lies close to the TAMS line for the most massiv e BSs and 
somewhat higher for TO stars (ISchaller et al.lll9^ . A similar 
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Table 2. Model binaries with periods shorter than 1 d at the age of 11 Gyr 


model 

Ml (Mo) 

M2(Mo) 

q 

F(d) log(L/Lo) 

logT, 

Remarks 

0.75-h0.4(1.6) 

0.458 

0.618 

0.741 

0.168 

-0.489 

3.723 

CB,MS 

0.75+0.5(1.5) 

0.526 

0.643 

0.818 

0.177 

-0.380 

3.731 

CB,MS 

0.75+0.7(1.5) 

0.580 

0.764 

0.759 

0.194 

-0.149 

3.757 

CB,MS 

0.8+0.4(1.8) 

0.441 

0.678 

0.650 

0.173 

-0.363 

3.737 

CB,MS 

0.8+0.5(1.7) 

0.544 

0.667 

0.816 

0.226 

-0.295 

3.737 

CB,MS 

0.8+0.6(1.6) 

0.547 

0.752 

0.727 

0.220 

-0.154 

3.754 

CB,MS 

0.8+0.7(1.6) 

0.535 

0.850 

0.629 

0.213 

0.017 

3.777 

CB,MS 

0.81+0.5(1.7) 

0.487 

0.732 

0.665 

0.201 

-0.223 

3.750 

CB,MS 

0.81+0.6(1.6) 

0.489 

0.818 

0.598 

0.198 

-0.057 

3.770 

CB,MS 

0.82+0.4(1.9) 

0.447 

0.690 

0.648 

0.183 

-0.327 

3.740 

CB,MS 

0.82+0.5(1.7) 

0.460 

0.767 

0.600 

0.166 

-0.166 

3.759 

CB,MS 

0.82+0.7(1.7) 

0.659 

0.742 

0.888 

0.283 

-0.085 

3.751 

CB,MS 

0.83+0.4(1.9) 

0.432 

0.713 

0.606 

0.145 

-0.284 

3.745 

CB,MS 

0.83+0.5(1.8) 

0.525 

0.710 

0.739 

0.222 

-0.217 

3.746 

CB,MS 

0.83+0.6(1.7) 

0.552 

0.771 

0.716 

0.240 

-0.101 

3.758 

CB,MS 

0.83+0.7(1.7) 

0.584 

0.827 

0.706 

0.256 

0.008 

3.769 

CB,MS 

0.84+0.4(2.0) 

0.442 

0.711 

0.622 

0.193 

-0.282 

3.745 

CB,MS 

0.84+0.5(1.8) 

0.451 

0.792 

0.569 

0.205 

-0.122 

3.765 

CB,MS 

0.84+0.6(1.7) 

0.481 

0.852 

0.565 

0.219 

0.011 

3.778 

CB,MS 

0.84+0.7(1.7) 

0.521 

0.898 

0.580 

0.229 

0.117 

3.788 

CB,MS 

0.85+0.5(1.9) 

0.542 

0.710 

0.763 

0.237 

-0.192 

3.746 

CB,MS 

0.85+0.6(1.8) 

0.611 

0.731 

0.836 

0.271 

-0.108 

3.750 

CB,MS 

0.86+0.4(2.1) 

0.415 

0.755 

0.550 

0.190 

-0.217 

3.755 

CB,MS 

0.86+0.5(1.9) 

0.452 

0.808 

0.559 

0.214 

-0.086 

3.769 

CB,MS 

0.86+0.6(1.8) 

0.507 

0.843 

0.601 

0.244 

0.010 

3.775 

CB,MS 

0.86+0.7(1.8) 

0.589 

0.847 

0.695 

0.275 

0.055 

3.772 

CB,MS 

0.865+0.4(2.2) 

0.474 

0.700 

0.677 

0.226 

-0.261 

3.743 

CB,MS 

0.865+0.5(1.9) 

0.444 

0.821 

0.541 

0.176 

-0.055 

3.772 

CB,MS 

0.865+0.5(2.0) 

0.556 

0.709 

0.784 

0.261 

-0.172 

3.746 

CB,MS 

0.865+0.6(1.8) 

0.478 

0.875 

0.546 

0.214 

0.063 

3.784 

CB,MS 

0.865+0.7(1.8) 

0.527 

0.913 

0.577 

0.260 

0.148 

3.790 

CB,MS 

0.865+0.8(1.9) 

0.662 

0.863 

0.767 

0.305 

0.117 

3.773 

CB,MS 

0.8675+0.4(2.2) 

0.378 

0.797 

0.474 

0.242 

-0.097 

3.767 

CB,MS 

0.8675+0.4(2.3) 

0.432 

0.744 

0.581 

0.292 

-0.123 

3.752 

CB,MS 

0.8675+0.5(2.0) 

0.420 

0.846 

0.496 

0.276 

0.025 

3.776 

CB,MS 

0.8675+0.5(2.1) 

0.508 

0.758 

0.670 

0.345 

-0.023 

3.754 

CB,MS 

0.8675+0.6(1.9) 

0.475 

0.879 

0.540 

0.321 

0.123 

3.780 

CB,MS 

0.8675+0.6(2.0) 

0.595 

0.759 

0.784 

0.412 

0.079 

3.754 

CB,MS 

0.8675+0.7(1.8) 

0.426 

1.012 

0.421 

0.242 

0.336 

3.818 

CB,BS 

0.8675+0.7(1.9) 

0.537 

0.903 

0.595 

0.372 

0.204 

3.781 

CB,MS 

0.8675+0.8(1.9) 

0.479 

1.046 

0.458 

0.331 

0.403 

3.821 

CB,BS 

0.8675+0.8(2.0) 

0.593 

0.934 

0.635 

0.489 

0.318 

3.780 

CB,MS 

0.88+0.5(2.0) 

0.191 

1.072 

0.178 

0.845 

0.543 

3.823 

NCB,BS 

0.9+0.4(2.4) 

0.189 

1.001 

0.189 

0.843 

0.436 

3.802 

NCB,MS 

0.9+0.6(1.9) 

0.152 

1.221 

0.124 

0.709 

0.836 

3.887 

NCB,BS 


red bou ndar y was adopted by lLeigh et al.l (l2007h and lLeigh et al.l 
(1201 lab . but lFiorentino et aLl(l2014h . for instance, used a vertical 
line in the HRD to delimit the BS region. 

Forty-five binaries (4.6% of all models) fulfill these con¬ 
ditions. Three of them have periods between 0.5 and 1 d, and 
the rest has periods shorter than 0.5 d. All CBs and NCBs are 
clustered near the TO point. However, most of them lie below 
TO, with a mere four falling in the region of BS. The situation 
is different for binaries with F > 1 d (not plotted). They all are 
of Algol-type with low mass ratios and lie in the BS region or 
redward of it. It is interesting to note that progenitors of CBs 
and NCBs have initial periods from a narrow interval: 33 out of 
45 had periods between 1.7-2.0 d and none had a period longer 
than 2.4 d. Their total initial masses are between 1.15-1.5 Mq 
, but they are lowered to 1.08-1.37 Mq at the age of 11 Gyr. 
Short periods dominate among the plotted binaries: apart from 
three stars with periods longer than 0.5 d (plus signs), there are 
only seven more with periods between 0.3-0.5 d and the rest. 


that is, 35 stars, has periods shorter than 0.3 d, with the record 
shortest equal to ~ 0.14 d. Similarly short periods dominate 
among field C Bs, although he re the shortest known period is 
close to 0.2 d (lRucinskill2007h . We can expect periods shorter 
than 0.2 d among low-metallicity CBs, however, because their 
components are more compact than their solar-metallicity coun¬ 
terparts. Binaries with periods shorter than 0.3 d have a very high 
mass ratio, between 0.5-0.8, whereas those with longer periods 
have lower ratios down to 0.1-0.2. This trend results from the 
dependence of the mass ratio past mass exchange on the differ¬ 
ence in core helium content between the components. For two 
chemically uniform ZA MS stars with un equal masses, the mass 
transfer stops at ^ = 1 iKuinetl (1194Ih . TTie difference in the 
core helium content between a more and less massive compo¬ 
nent increases in the course of evolution, so that the later mass 
transfer takes place, the lower the mass ratio. Observations of 
field CBs with periods shorter than 0.3 d indicate that most of 
them belong to least massive CBs with both components still in 
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Table 3. BS mergers formed between 10 and 11 Gyr, evolved to 
the age of 11 Gyr 


model rmera(years) M{Mq) log(L/LQ) log 

0.75+0.4(1.5) 1 


0.8+0.4(1.7) 1 

0.8+0.5(1.6) 1 

0.8+0.6(1.5) 1 

0.8+0.7(1.5) 1 

0.81+0.4(1.7) 1 

0.81+0.4(1.8) 1 

0.81+0.5(1.6) 1 

0.81+0.6(1.5) 1 

0.82+0.4(1.8) 1 

0.82+0.5(1.6) 1 

0.82+0.6(1.6) 1 

0.83+0.4(1.8) 1 

0.83+0.5(1.7) 1 

0.83+0.6(1.6) 1 

0.84+0.4(1.9) 1 

0.84+0.5(1.7) 1 

0.84+0.6(1.6) 1 

0.85+0.4(1.9) 1 

0.85+0.4(2.0) 1 

0.85+0.5(1.7) 1 

0.85+0.5(1.8) 1 

0.85+0.6(1.7) 1 

0.85+0.7(1.7) 1 

0.86+0.4(2.0) 1 

0.86+0.5(1.8) 1 

0.86+0.6(1.7) 1 

0.86+0.7(1.7) 1 

0.865+0.4(2.0) 1 

0.865+0.4(2.1) 1 

0.865+0.5(1.8) 1 

0.865+0.6(1.7) 1 

0.865+0.8(1.8) 1 

0.8675+0.4(2.0) 1 

0.8675+0.4(2.1) 1 

0.8675+0.5(1.8) 1 

0.8675+0.5(1.9) 1 

0.8675+0.6(1.8) 1 

0.87+0.4(2.0) 1 

0.87+0.5(1.9) 1 

0.88+0.4(2.1) 1 

0.88+0.5(1.9) 1 

0 . 88 + 0 . 6 ( 1 . 8 ) 1 

0.9+0.4(2.2) 1 


.037E+10 

0.981 

.040E+10 

1.024 

.075E+10 

1.114 

.042E+10 

1.204 

.OllE+10 

1.294 

.005E+10 

1.033 

.092E+10 

1.028 

.044E+10 

1.122 

.014E+10 

1.216 

.057E+10 

1.039 

.014E+10 

1.134 

.086E+10 

1.217 

.018E+10 

1.051 

.072E+10 

1.138 

.055E+10 

1.229 

.059E+10 

1.056 

.038E+10 

1.149 

.024E+10 

1.241 

.019E+10 

1.068 

.099E+10 

1.061 

.002E+10 

1.160 

.087E+10 

1.153 

.085E+10 

1.242 

.085E+10 

1.330 

.049E+10 

1.074 

.046E+10 

1.165 

.045E+10 

1.254 

.053E+10 

1.340 

.023E+10 

1.080 

.099E+10 

1.074 

.021E+10 

1.172 

.022E+10 

1.262 

.064E+10 

1.429 

.015E+10 

1.083 

.090E+10 

1.077 

.OllE+10 

1.173 

.086E+10 

1.168 

.093E+10 

1.254 

.005E+10 

1.086 

.080E+10 

1.168 

.031E+10 

1.092 

.033E+10 

1.181 

.041E+10 

1.271 

.OllE+10 

1.110 


0.331 

3.841 

0.467 

3.857 

0.657 

3.888 

0.971 

3.912 

1.072 

3.906 

0.529 

3.862 

0.439 

3.856 

0.728 

3.891 

1.141 

3.793 

0.505 

3.862 

0.846 

3.894 

0.856 

3.910 

0.572 

3.869 

0.734 

3.895 

0.995 

3.913 

0.547 

3.869 

0.830 

3.898 

1.251 

3.748 

0.625 

3.875 

0.528 

3.868 

1.003 

3.863 

0.749 

3.899 

0.911 

3.912 

1.048 

3.921 

0.615 

3.876 

0.867 

3.905 

1.128 

3.831 

1.341 

3.733 

0.672 

3.879 

0.566 

3.874 

0.987 

3.885 

1.133 

3.829 

1.802 

3.716 

0.728 

3.881 

0.637 

3.878 

1.017 

3.773 

0.827 

3.897 

0.983 

3.916 

0.772 

3.880 

0.897 

3.908 

0.822 

3.878 

1.029 

3.760 

1.275 

3.739 

0.986 

3.808 


the early phases of the MS life (IStenieri & Gazeaj|2012h . Their 
mass ratios clus ter around 0.5, as opposed to longer period CBs 
(lRucinskill2010l) . Most CBs with periods shorter than 0.3 d stay 
in contact for less than 1 Gyr, and then coalescence occurs as a 
consequence of the outer Roche lobe overflow. The longer pe- 
riod binaries liv e longer as CBs or NCBs (up to 2 Gyr, see also 
iTian et aL l|2006b before merging as a result of the Darwin insta- 
bilitv (iRasiol 199^ . 

In addition to binaries, 44 mergers (4.5% of all models), 
formed between 10 and 11 Gyr and evolved to the age of 11 
Gyr, are listed in Table [3] and plotted in Fig. as open cir¬ 
cles. The columns give the name of the model, age of merging, 
merger mass, luminosity, and effective temperature. In addition 
to 34 stars lying in the BS region, 10 subgiants or giants are 
also shown for comparison. These are mergers formed shortly 
after 10 Gyr, which already left the MS and now approach the 
red giant branch. Because the orbital AM of a binary just be¬ 
fore merging is far too high to be accommodated in a spin AM 



Fig. 6. HR diagram with model CBs, NCBs, and mergers at 
the age of 11 Gyr, together with the corresponding low-metal 
isochrone. Asterisks denote binaries with an orbital period 
shorter than 0.5 d, plus signs binaries with periods between 0.5 
and 1 d, and open circles denote mergers formed past 10 Gyr 
and evolved to the age of 11 Gyr. The ZAMS is marked with a 
dotted line and the adopted red boundary of the BS region with 
long dashes. 


of a single star, most of it must be expelled together with some 
amount of mass in t he form of a ring or dis k surrounding the 
freshly formed star dRasio & ShaDiroill995h . We adopted 0.1 
Mq for th e expelled mass, based on an indirect argument pre¬ 
sented by IStenieh & Kiragal (l2013bh . This ma ss ejected to the 
distance of 10-20 Rn can absorb all excess AM (ISills et al.l200ll: 
IStenieh & Kiragall2013bh . The resulting merger mass is equal to 
the sum of the component masses minus 0.1 Mq. Mergers are as¬ 
sumed to evolve like single stars with the same degree of hydro¬ 
gen depletion (L_e^ no mixing of the merger matter is assumed, 
see ISills et al.ll2009l) . 

Disregarding subgiants, mergers form a well-defined se¬ 
quence that lies close to the ZAMS for the lowest masses and 
increasingly deviates from it for the increasing mass. The most 
massive merger with a mass of 1.43 Mq also has the highest lu¬ 
minosity and already climbs the red giant branch (Fig. IDl. All 
other mergers have masses between 0.98-1.33 Mq. 

How does the situation change with age? Figure [T] shows 
similar data as Fig. 0 but 1 Gyr later. Here, again, all binaries 
with periods shorter than 1 d at the age of 12 Gyr are plotted 
with the same coding as in Fig. |6] Most binaries with periods 
shorter than 0.5 d are new because almost all those plotted in 
Fig. |6]did not survive till 12 Gyr and merged in the meantime. 
But seven binaries from Fig.|6]survived for more than 1 Gyr and 
are also plotted in Fig. [T] Mergers formed between 11-12 Gyr are 
shown as open circles. They are supplemented with the mergers 
from the previous figure, which did not yet climb the red giant 
branch. Those mergers are shown as filled circles. The numerical 
data are given in Tables |4] and |5] 
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Table 4. Model binaries with periods shorter than 1 d at the age of 12 Gyr 


model 

Ml (Mo) 

M2(Mo) 

q 

F(d) 

log(L/Lo) 

logT, 

Remarks 

0.7+0.4(1.5) 

0.445 

0.583 

0.763 

0.142 

-0.528 

3.720 

CB,MS 

0.75-r0.5(1.6) 

0.541 

0.619 

0.874 

0.199 

-0.398 

3.728 

CB,MS 

0.75-r0.7(1.6) 

0.603 

0.731 

0.825 

0.223 

-0.194 

3.749 

CB,MS 

0.8+0.4(1.9) 

0.426 

0.686 

0.621 

0.167 

-0.348 

3.739 

CB,MS 

0.8-r0.5(1.8) 

0.545 

0.657 

0.830 

0.229 

-0.287 

3.736 

CB,MS 

0.8-r0.6(1.7) 

0.552 

0.738 

0.748 

0.229 

-0.161 

3.751 

CB,MS 

0.8-r0.7(1.7) 

0.533 

0.842 

0.633 

0.231 

0.002 

3.775 

CB,MS 

0.81+0.5(1.8) 

0.468 

0.743 

0.630 

0.206 

-0.203 

3.752 

CB,MS 

0.81+0.6(1.7) 

0.471 

0.827 

0.570 

0.215 

-0.044 

3.773 

CB,MS 

0.81+0.7(1.7) 

0.511 

0.872 

0.586 

0.187 

0.070 

3.782 

CB,MS 

0.82+0.4(2.0) 

0.410 

0.718 

0.571 

0.183 

-0.290 

3.746 

CB,MS 

0.82+0.5(1.8) 

0.439 

0.780 

0.563 

0.162 

-0.141 

3.762 

CB,MS 

0.82+0.6(1.8) 

0.602 

0.705 

0.854 

0.280 

-0.155 

3.746 

CB,MS 

0.82+0.7(1.8) 

0.634 

0.759 

0.835 

0.295 

-0.066 

3.754 

CB,MS 

0.83+0.4(2.1) 

0.458 

0.679 

0.675 

0.208 

-0.312 

3.738 

CB,MS 

0.83+0.5(1.9) 

0.486 

0.742 

0.655 

0.223 

-0.183 

3.752 

CB,MS 

0.83+0.6(1.8) 

0.526 

0.788 

0.668 

0.248 

-0.071 

3.762 

CB,MS 

0.83+0.7(1.8) 

0.565 

0.836 

0.676 

0.264 

0.026 

3.771 

CB,MS 

0.84+0.4(2.1) 

0.398 

0.747 

0.533 

0.183 

-0.236 

3.753 

CB,MS 

0.84+0.5(1.9) 

0.424 

0.811 

0.523 

0.198 

-0.086 

3.770 

CB,MS 

0.84+0.6(1.8) 

0.457 

0.866 

0.528 

0.220 

0.041 

3.782 

CB,MS 

0.84+0.7(1.8) 

0.498 

0.911 

0.547 

0.236 

0.138 

3.791 

CB,MS 

0.85+0.4(2.2) 

0.435 

0.717 

0.607 

0.209 

-0.252 

3.746 

CB,MS 

0.85+0.5(2.0) 

0.487 

0.756 

0.644 

0.239 

-0.141 

3.755 

CB,MS 

0.85+0.6(1.9) 

0.565 

0.766 

0.738 

0.273 

-0.067 

3.756 

CB,MS 

0.85+0.7(1.9) 

0.650 

0.769 

0.845 

0.312 

-0.007 

3.756 

CB,MS 

0.86+0.4(2.2) 

0.391 

0.772 

0.506 

0.166 

-0.178 

3.760 

CB,MS 

0.86+0.4(2.3) 

0.459 

0.704 

0.652 

0.228 

-0.248 

3.744 

CB,MS 

0.86+0.5(2.0) 

0.420 

0.834 

0.504 

0.200 

-0.035 

3.776 

CB,MS 

0.86+0.5(2.1) 

0.537 

0.714 

0.752 

0.277 

-0.159 

3.747 

CB,MS 

0.86+0.6(1.9) 

0.466 

0.873 

0.534 

0.248 

0.064 

3.783 

CB,MS 

0.86+0.7(1.9) 

0.543 

0.883 

0.615 

0.279 

0.109 

3.781 

CB,MS 

0.865+0.4(2.3) 

0.413 

0.753 

0.548 

0.214 

-0.197 

3.754 

CB,MS 

0.865+0.5(2.1) 

0.485 

0.771 

0.629 

0.256 

-0.103 

3.758 

CB,MS 

0.865+0.6(1.9) 

0.452 

0.891 

0.507 

0.196 

0.100 

3.789 

CB,MS 

0.865+0.6(2.0) 

0.587 

0.755 

0.777 

0.314 

-0.030 

3.754 

CB,MS 

0.865+0.7(1.9) 

0.485 

0.944 

0.514 

0.267 

0.204 

3.799 

CB,MS 

0.865+0.8(2.0) 

0.574 

0.939 

0.611 

0.311 

0.221 

3.791 

CB,MS 

0.8675+0.4(2.3) 

0.215 

0.952 

0.226 

0.388 

0.229 

3.804 

CB,MS 

0.8675+0.4(2.4) 

0.233 

0.931 

0.250 

0.553 

0.244 

3.790 

NCB,MS 

0.8675+0.4(2.5) 

0.275 

0.887 

0.310 

0.705 

0.262 

3.775 

NCB,MS 

0.8675+0.4(2.6) 

0.283 

0.877 

0.323 

0.955 

0.352 

3.768 

NCB,MS 

0.8675+0.5(2.1) 

0.221 

1.032 

0.214 

0.583 

0.418 

3.819 

NCB,BS 

0.8675+0.5(2.2) 

0.279 

0.972 

0.287 

0.711 

0.377 

3.792 

NCB,MS 

0.8675+0.6(1.9) 

0.179 

1.161 

0.154 

0.579 

0.668 

3.872 

NCB,BS 

0.8675+0.6(2.0) 

0.256 

1.080 

0.237 

0.790 

0.552 

3.819 

NCB,BS 

0.8675+0.7(1.9) 

0.224 

1.199 

0.187 

0.655 

0.750 

3.873 

NCB,BS 

0.8675+0.7(2.0) 

0.296 

1.126 

0.263 

0.924 

0.645 

3.821 

NCB,BS 

0.8675+0.8(1.9) 

0.148 

1.363 

0.109 

0.630 

1.078 

3.920 

NCB,BS 

0.8675+0.8(2.0) 

0.275 

1.229 

0.224 

0.940 

0.816 

3.858 

NCB,BS 


The CBs and NCBs present at the age of 12 Gyr have similar 
properties as those existing 1 Gyr earlier. Almost all short-period 
binaries lie below or very close to the TO point. Out of 39 bina¬ 
ries (4% of all modeled), only one lies slightly above TO, but 
out of 11 long-period binaries 7 lie in the region of BS (two are 
hidden in a dense clump of filled circles). Similarly as in the pre¬ 
vious case, most of the initial orbital periods come from a narrow 
interval of 1.8-2.1 d that is shifted longward by 0.1 d compared 
to the binaries plotted in Fig. |6] The total initial masses lie be¬ 
tween 1.1-1.67 Mq , but they are lowered to 1.03-1.5 Mq at 12 
Gyr. 

Of the mergers formed between 11-12 Gyr (open circles in 
Fig. |7]i, 30 (3% of the total number of CCBMs) lie in the BS 


region and 11 (1.1% ) are subgiants. An additional 16 mergers 
(1.6%) were formed before 11 Gyr (filled circles). Twelve of 
these still lie in the BS region, which indicates that a significant 
fraction of BS spends more than 1 Gyr in this region. Three of 
them still remain in the BS region even at the age of 13 Gyr, that 
is, almost 3 Gyr after formation (see below). 

The situation after another 1 Gyr is displayed in Fig.[8]with 
the same coding as in Figs.|6]l2] The numerical data of the plotted 
models are given in Tables|2and|7] We have here 46 short-period 
(4.7%) and 6 long-period (0.6%) CBs and NCBs. Only four bi¬ 
naries (two short-period and two long-period binaries) lie above 
the 13 Gyr isochrone in the region of BS, and all others lie below 
the TO point. Similarly as in the earlier age, the initial periods of 
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Fig. 7. HR diagram with model CBs, NCBs, and mergers at 
the age of 12 Gyr, together with the corresponding low-metal 
isochrone. Asterisks denote binaries with an orbital period 
shorter than 0.5 d, plus signs binaries with periods between 0.5 
and 1 d, filled circles denote mergers formed between 10-11 Gyr, 
and open circles those formed between 11-12 Gyr. All mergers 
are evolved to the age of 12 Gyr. The ZAMS is marked with a 
dotted line and the adopted red boundary of the BS region with 
long dashes. 


most of the binaries come from a narrow range 1.9-2.2 d, shifted 
again longward by 0.1 d, compared to the age of 12 Gyr. The 
initial total masses of the binaries are between 1.2-1.67 Mq , but 
the present masses are lowered to 1.0-1.5 Mq. 

Thirty-six mergers (3.7%), formed past 12 Gyr, are plot¬ 
ted in Fig. [8] as open circles. They are supplemented with 24 
mergers (2.5%) formed between 11 and 12 Gyr and 3 more, 
formed shortly after 10 Gyr. Fifteen mergers fall in the subgiant 
region, and the rest lies in the BS region, among them the three 
low-massbinaries mentioned above. Merger masses do not differ 
from the masses of the mergers plotted in the previous figures: 
they fall into the interval of 0.93-1.32 Mq. 

4. Comparison of model calculations with 
observations 

4.1. Contact and near-contact binaries 

Recently, results of a deep search for photometric variables in 
a few GCs have been publ ished dKaluznv & Thompsonll200^ 
iKaluznv et alll2010l l2013all^ . In all cases, no CBs were de¬ 
tected among unevolved cluster members, although a prominent 
binary se quence is visible abo ve the cluster MS on the CMD 
(see e. g.. IKaluznv et al.ll2013al) . The authors concluded that CBs 
begin only to appear about 1 mag below the TO. Our results 
are fully consistent with this conclusion. Assuming that TO lies 
approximately at logL = 0 in Figs. IMSl the faintest CBs lie 
about 0.4-0.5 in logL below TO, which corresponds to 1-1.2 
mag. All model binaries with still lower luminosity are detached. 



Fig. 8. HR diagram with model CBs, NCBs, and mergers at 
the age of 13 Gyr, together with the corresponding low-metal 
isochrone. Asterisks denote binaries with an orbital period 
shorter than 0.5 d, plus signs binaries with periods between 0.5 
and 1 d, filled circles denote mergers formed between 10-12 Gyr, 
and open circles those formed between 12-13 Gyr. All mergers 
are evolved to the age of 13 Gyr. The ZAMS is marked with a 
dotted line and the adopted red boundary of the BS region with 
long dashes. 


The complete lack of faint, low-mass model CBs results from 
the adopted value of 1.5 d for the initial orbital cut-off period. 
Lowering this value (e. g., to 1 d or less) leads to a contact con¬ 
figuration reached within the age of a GC by binaries with essen¬ 
tially unevolved, low-mass components. Their absence among 
the observed CBs is direct evidence that the cut -off period is no t 
shorter than 1.5 d. It may even be close to 2 d (IStepiehll201 lal) . 
Adopting the latter value would result in a disappearance of sev¬ 
eral of the faintest CBs in Figs. |6]l8] (and some of the plotted 
mergers). 

On the whole, of 975 CCBMs, only about 50 (~ 5%) 
are in a contact phase at a given age (more accurately, there 
are 45, 50, and 52 at the age of 11, 12, and 13 Gyr, respec¬ 
tively). Approximately 10% of these CBs fall in the BS region. 
Qualit atively similar proportions were obtained by iTian et aP 
(|2006|) . They modeled over one million binaries, and only less 
than 2000 were in the mass-exchange phase at the age of 4 Gyr, 
corresponding to the age of M67. About 11% of these fell in 
the BS region. A more quantitative comparison between their 
and our results cannot be made because of distinctly different 
assumptions adopted. 

A lower CB frequency in GC s than in the field has already 
been noticed by iRucinsO (l2000h . It results most likely from a 
yery low binary frequency of ~ 0.1 in GCs (but with a signif¬ 
icant scatter among diffe rent clusters), compared to 0.5 in the 
field dMilone et alJl2012li . together with a low percentage of bi¬ 
naries in the contact phase at a giyen age, as shown aboye. Any 
statistical comparison of models with obseryations is therefore 
extremely uncertain because we haye so few obseryational data. 
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Table 6. Model binaries with periods shorter than 1 d at the age of 13 Gyr 


model 

Ml (Mo) 

M2(Mo) 

q 

F(d) log(L/Lo) 

logT, 

Remarks 

0.7+0.5(1.5) 

0.480 

0.633 

0.758 

0.142 

-0.418 

3.729 

CB,MS 

0.7-r0.6(1.5) 

0.518 

0.680 

0.762 

0.178 

-0.326 

3.738 

CB,MS 

0.75-r0.6(1.6) 

0.509 

0.732 

0.695 

0.185 

-0.230 

3.750 

CB,MS 

0.8-r0.4(2.0) 

0.408 

0.699 

0.584 

0.167 

-0.328 

3.742 

CB,MS 

0.8-r0.5(1.9) 

0.528 

0.667 

0.792 

0.236 

-0.276 

3.738 

CB,MS 

0.8-r0.6(1.8) 

0.565 

0.716 

0.789 

0.241 

-0.171 

3.747 

CB,MS 

0.8-r0.7(1.8) 

0.551 

0.816 

0.675 

0.245 

-0.026 

3.768 

CB,MS 

0.81+0.4(2.1) 

0.446 

0.668 

0.668 

0.199 

-0.338 

3.736 

CB,MS 

0.81+0.5(1.9) 

0.450 

0.754 

0.597 

0.210 

-0.186 

3.755 

CB,MS 

0.81+0.6(1.8) 

0.468 

0.822 

0.569 

0.225 

-0.047 

3.771 

CB,MS 

0.81+0.7(1.8) 

0.493 

0.881 

0.560 

0.210 

0.080 

3.785 

CB,MS 

0.82+0.4(2.1) 

0.386 

0.737 

0.524 

0.179 

-0.256 

3.751 

CB,MS 

0.82+0.5(1.9) 

0.410 

0.802 

0.511 

0.170 

-0.101 

3.769 

CB,MS 

0.82+0.5(2.0) 

0.537 

0.675 

0.796 

0.253 

-0.220 

3.740 

CB,MS 

0.82+0.6(1.9) 

0.589 

0.709 

0.831 

0.287 

-0.132 

3.746 

CB,MS 

0.82+0.7(1.9) 

0.636 

0.747 

0.851 

0.307 

-0.048 

3.752 

CB,MS 

0.83+0.4(2.2) 

0.413 

0.717 

0.576 

0.204 

-0.266 

3.746 

CB,MS 

0.83+0.5(2.0) 

0.448 

0.771 

0.581 

0.227 

-0.141 

3.759 

CB,MS 

0.83+0.6(1.9) 

0.500 

0.806 

0.620 

0.255 

-0.043 

3.766 

CB,MS 

0.83+0.7(1.9) 

0.541 

0.849 

0.637 

0.277 

0.051 

3.774 

CB,MS 

0.84+0.4(2.2) 

0.375 

0.763 

0.491 

0.163 

-0.199 

3.758 

CB,MS 

0.84+0.4(2.3) 

0.449 

0.690 

0.651 

0.228 

-0.264 

3.741 

CB,MS 

0.84+0.5(2.0) 

0.398 

0.830 

0.480 

0.192 

-0.045 

3.776 

CB,MS 

0.84+0.5(2.1) 

0.521 

0.707 

0.737 

0.263 

-0.168 

3.746 

CB,MS 

0.84+0.6(1.9) 

0.429 

0.884 

0.485 

0.222 

0.077 

3.787 

CB,MS 

0.84+0.6(2.0) 

0.611 

0.704 

0.868 

0.324 

-0.084 

3.746 

CB,MS 

0.84+0.7(1.9) 

0.465 

0.934 

0.498 

0.249 

0.182 

3.798 

CB,MS 

0.85+0.4(2.3) 

0.383 

0.764 

0.501 

0.204 

-0.191 

3.758 

CB,MS 

0.85+0.4(2.4) 

0.461 

0.685 

0.673 

0.265 

-0.246 

3.741 

CB,MS 

0.85+0.5(2.1) 

0.438 

0.795 

0.551 

0.240 

-0.085 

3.765 

CB,MS 

0.85+0.6(2.0) 

0.514 

0.809 

0.635 

0.281 

-0.011 

3.766 

CB,MS 

0.85+0.7(2.0) 

0.600 

0.808 

0.743 

0.322 

0.047 

3.763 

CB,MS 

0.86+0.4(2.4) 

0.402 

0.753 

0.534 

0.221 

-0.190 

3.755 

CB,MS 

0.86+0.5(2.2) 

0.486 

0.756 

0.643 

0.270 

-0.109 

3.755 

CB,MS 

0.86+0.6(2.0) 

0.422 

0.908 

0.465 

0.250 

0.127 

3.793 

CB,MS 

0.86+0.6(2.1) 

0.612 

0.717 

0.854 

0.333 

-0.038 

3.748 

CB,MS 

0.86+0.7(2.0) 

0.498 

0.916 

0.544 

0.289 

0.168 

3.790 

CB,MS 

0.865+0.4(2.4) 

0.346 

0.811 

0.427 

0.220 

-0.092 

3.771 

CB,MS 

0.865+0.4(2.5) 

0.405 

0.753 

0.538 

0.274 

-0.137 

3.754 

CB,MS 

0.865+0.4(2.6) 

0.398 

0.760 

0.524 

0.414 

-0.041 

3.754 

CB,MS 

0.865+0.5(2.2) 

0.403 

0.844 

0.477 

0.267 

0.012 

3.776 

CB,MS 

0.865+0.5(2.3) 

0.496 

0.750 

0.661 

0.342 

-0.039 

3.753 

CB,MS 

0.865+0.5(2.4) 

0.523 

0.725 

0.721 

0.524 

0.084 

3.750 

NCB,MS 

0.865+0.6(2.1) 

0.478 

0.857 

0.558 

0.328 

0.091 

3.775 

CB,MS 

0.865+0.6(2.2) 

0.508 

0.827 

0.614 

0.516 

0.158 

3.763 

NCB,MS 

0.865+0.7(2.0) 

0.411 

1.007 

0.408 

0.281 

0.322 

3.817 

CB,BS 

0.865+0.7(2.1) 

0.501 

0.916 

0.547 

0.451 

0.254 

3.781 

CB,MS 

0.865+0.8(2.1) 

0.470 

1.031 

0.456 

0.330 

0.378 

3.818 

CB,BS 

0.865+0.7(2.2) 

0.627 

0.792 

0.792 

0.617 

0.253 

3.757 

NCB,MS 

0.865+0.8(2.2) 

0.527 

0.975 

0.541 

0.584 

0.389 

3.786 

NCB,MS 

0.8675+0.4(2.4) 

0.125 

1.023 

0.122 

0.806 

0.443 

3.818 

NCB,BS 

0.8675+0.5(2.1) 

0.104 

1.133 

0.092 

0.875 

0.691 

3.860 

NCB,BS 


Nevertheless, we can compare the predicted to observed frac¬ 
tion of BSs among all CBs in a GC. We used data for two clus¬ 
ters observed by Kaluzny and his collaborators w ith the highest 
number of CBs (except for w Centauri). In M4, iKaluznv et alJ 
(l2Q13a) detected nine CBs with one BS among them, and in 47 
Tuc. Kaluzny et al](l2013bh identified 15 CBs orNCBs, of which 
six are BS. The resulting fractions are 0.1 and 0.4, respectively. 
We di d not take into accou nt w Cen with the richest population of 
CBs (IKaluznv et al.ll2004l) because it is highly atypical and the 
accurate number of member CBs is not well known, although 
the approximate data indicate a similarly high ratio as in case of 


47 Tuc. The predicted fractions range from 0.08 to 0.2 for the 
three considered cluster ages. As we see, they are close to the 
lower observed value but are at odds with the higher value. This 
may indicate some deficiencies of the CCBM, for instance, too 
short initial cut-off periods, a too high AML rate at short peri¬ 
ods, or a too low mass-transfer rate. With the lower AML rate 
and/or higher mass transfer rate, a binary stays longer in con¬ 
tact and can reach a lower mass ratio before merging. Many of 
the short-period CBs have a rather high mass ratio of about 0.7- 
0.8 at the time of merging, whereas field CBs with P % 0.3 d 
center around a value of 0.5 (lRucinskill2010t IStenieh & GazeasI 
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Table 5. BS mergers formed between 10 and 12 Gyr, evolved to 
the age of 12 Gyr 


model rmera(years) M{Mq) log(L/LQ) log 

0.75+0.4(1.5) 1 


0.8+0.4(1.7) 1 

0.8+0.5(1.6) 1 

0.81+0.4(1.7) 1 

0.81+0.4(1.8) 1 

0.82+0.4(1.8) 1 

0.83+0.4(1.8) 1 

0.84+0.4(1.9) 1 

0.85+0.5(1.8) 1 

0.86+0.4(2.0) 1 

0.865+0.4(2.0) 1 

0.865+0.4(2.1) 1 

0.8675+0.4(2.0) 1 

0.8675+0.4(2.1) 1 

0.87+0.4(2.0) 1 

0.87+0.5(1.9) 1 

0.75+0.4(1.6) 1 

0.75+0.5(1.5) 1 

0.75+0.6(1.5) 1 

0.75+0.7(1.5) 1 

0.8+0.4(1.8) 1 

0.8+0.5(1.7) 1 

0 . 8 + 0 . 6 ( 1 . 6 ) 1 

0.81+0.4(1.9) 1 

0.81+0.5(1.7) 1 

0.81+0.6(1.6) 1 

0.82+0.4(1.9) 1 

0.82+0.5(1.7) 1 

0.82+0.6(1.7) 1 

0.82+0.7(1.7) 1 

0.83+0.4(1.9) 1 

0.83+0.4(2.0) 1 

0.83+0.5(1.8) 1 

0.83+0.6(1.7) 1 

0.83+0.7(1.7) 1 

0.84+0.4(2.0) 1 

0.84+0.5(1.8) 1 

0.84+0.6(1.7) 1 

0.85+0.4(2.1) 1 

0.85+0.5(1.9) 1 

0.85+0.6(1.8) 1 

0.85+0.7(1.8) 1 

0.86+0.4(2.1) 1 

0.86+0.5(1.9) 1 

0 . 86 + 0 . 6 ( 1 . 8 ) 1 

0.86+0.7(1.8) 1 

0.865+0.4(2.2) 1 

0.865+0.5(1.9) 1 

0.865+0.5(2.0) 1 

0.8675+0.4(2.2) 1 

0.8675+0.5(2.0) 1 

0.87+0.4(2.1) 1 

0.87+0.4(2.2) 1 

0.87+0.6(1.8) 1 

0.88+0.4(2.2) 1 

0.89+0.4(2.2) 1 

0.89+0.4(2.3) 1 


.037E+10 

0.981 

.040E+10 

1.024 

.075E+10 

1.114 

.005E+10 

1.033 

.092E+10 

1.028 

.057E+10 

1.039 

.018E+10 

1.051 

.059E+10 

1.056 

.087E+10 

1.153 

.049E+10 

1.074 

.023E+10 

1.080 

.099E+10 

1.074 

.015E+10 

1.083 

.090E+10 

1.077 

.005E+10 

1.086 

.080E+10 

1.168 

.146E+10 

0.974 

.143E+10 

1.065 

.194E+10 

1.148 

.126E+10 

1.242 

.133E+10 

1.016 

.178E+10 

1.105 

.149E+10 

1.195 

.183E+10 

1.021 

.143E+10 

1.116 

.117E+10 

1.206 

.141E+10 

1.034 

.108E+10 

1.127 

.188E+10 

1.208 

.173E+10 

1.294 

.103E+10 

1.045 

.187E+10 

1.038 

.167E+10 

1.130 

.154E+10 

1.219 

.144E+10 

1.307 

.144E+10 

1.050 

.129E+10 

1.141 

.120E+10 

1.230 

.180E+10 

1.056 

.177E+10 

1.146 

.178E+10 

1.234 

.183E+10 

1.319 

.127E+10 

1.068 

.130E+10 

1.158 

.137E+10 

1.246 

.149E+10 

1.331 

.178E+10 

1.069 

.105E+10 

1.164 

.189E+10 

1.157 

.158E+10 

1.070 

.158E+10 

1.161 

.148E+10 

1.068 

.179E+10 

1.055 

.162E+10 

1.238 

.151E+10 

1.069 

.123E+10 

1.088 

.144E+10 

1.076 


0.418 3.846 

0.584 3.862 

0.864 3.881 

0.678 3.863 

0.546 3.862 

0.635 3.866 

0.743 3.867 

0.694 3.871 

1.054 3.781 

0.809 3.869 

0.940 3.817 

0.715 3.878 

1.468 3.721 

0.858 3.860 

0.965 3.845 

1.054 3.819 

0.327 3.840 

0.538 3.870 

0.683 3.898 

1.085 3.893 

0.476 3.856 

0.646 3.885 

0.945 3.913 

0.448 3.855 

0.733 3.889 

1.117 3.791 

0.521 3.862 

0.859 3.888 

0.845 3.907 

1.053 3.923 

0.595 3.868 

0.498 3.861 

0.743 3.893 

1.017 3.911 

1.461 3.726 

0.573 3.868 

0.873 3.894 

1.436 3.724 

0.559 3.869 

0.772 3.896 

0.965 3.918 

1.049 3.919 

0.651 3.875 

0.936 3.891 

1.206 3.754 

1.430 3.727 

0.600 3.875 

1.006 3.783 

0.784 3.898 

0.723 3.875 

0.966 3.886 

1.369 3.723 

0.910 3.754 

0.890 3.755 

1.152 3.739 

1.036 3.746 

1.110 3.741 


l2012h . The lower mass ratio means more mass transferred to the 
gainer, hence its higher mass and higher position on the HRD, 
that is, a higher probability for entering the BS region. There 
may also still be another explanation of the discrepancy: strong 
fluctuations of this fraction among different clusters suggest a 
nonuniform BSs formation rate in some clusters, with individ¬ 


ual bursts occuri'ing in the recent past (see below). Regardless of 
the reason, it is apparent that profound differences occur among 
GCs, which makes a comparison of theoretical predictions with 
individual clusters uncertain. 

A direct comparison of CCBMs to the observed param¬ 
eters of individual CBs or NCBs is hampered because only 
very few high-quality data are available for these GC mem¬ 
bers. To our knowledge, V60, an Algol-type BS in M55, is the 
only NCB with pararn eters based on light and velocity curves 
dRozvczka et alJl2013h . Generally, to hnd a progenitor for a 
given binary, a hne-tuning of the initial parameters is needed. 
Here we do not attempt it, but we can still hnd a model that re¬ 
produces V60 quite well. This is the most massive CCBM, that 
is, 0.89+0.8(1.9). Its parameters at the age of about 11 Gyr are 
compared with the observed values for V60 in Table [8] The po¬ 
sitions of the two binaries are also shown in the HRD (Fig. |9]l, 
where a diamond corresponds to V60 and the plus sign to the 
model. Another model, 0.9+0.8(1.9), evolves in about 11 Gyr to 
a similar conhguration with a period of 1.18 d and component 
masses of 0.24 and 1.30 Mq, which is also not far from V60, but 
its luminosity, logL = 0.998, is somewhat too high. 

As we see, the model parameters agree with the observed 
values within about 2cr, even though the model is not sufficiently 
massive to better reproduce the observations. The total mass of 
V60 is equal to 1.59 Mq, compared to only 1.54 Mq for the 
model. A slight mass increase, together with a correct period and 
metallicity should give a better agreement. When the initial pa¬ 
rameter values of CCBM were specihed, the following evolution 
was fully determined - we did not introduce any free parameters 
that would, for instance, describe additional mass and/or AM 
loss. 

The mass and radius of the less massive component of V60 
indicate that it is an evolved star with a helium core of ~ 0. 15Mq 
, whereas its companion is a ZAMS star. The binary is thus a typ¬ 
ical short-period Algol. It is commonly accepted that such a con¬ 
hguration results from a rapid mass transfer following a RLOF 
by an evolved primary. In effect, a binary is formed showing 
the well-known Algol paradox. No other way of producing an 
Algol-type binary seems to be viable, although many details of 
evolution preceding RLOF, starting from the origin of the ini¬ 
tial binary, possible dynamical interactions with other stars, or 
amount of mass and AM lost during the consecutive evolution¬ 
ary phases, are still unknown and need to be assumed for the cal¬ 
culated model. Even if all these details are hxed, as in case of the 
CCBM code, the results are equivocal. Different combinations of 
the initial parameters may res ult in the same hnal values, except 
that they have different age (ISteniehl 1201 Ibl : iPietrzvhski et al.l 
l2012h . If the age of the modeled binary is not known, ambiguity 
remains. For a binary member of a GC the age is known, which 
additionally constrains a possible initial model. 

Stellar parameters have also been obtained for a few more 
CBs and NCBs in GCs, based solely on the photo metric observa 
tions These are NH 19, 30 and 31 i n NGC 5466 JKa 


NIL 5 and V239 i n co Cen (iHelt et aljfl^ 


2012h . VI in NGC 6397 (iRubenstein & Bailvnlll996h . V228 in 


hath et al 


Li & Qian 


47 T uc (iKaluznv et al.ll2007t rSarnall2008h . and V47 and V53 in 
M4 (iLiu et al.ll20riir in most cases, however, the resulting bi¬ 
nary parameters are very uncertain, so any hne-tuning of the 
initial parameters to reproduce the obtained results must await 
better observations. Still, if the data are acceptable, we can al¬ 
ways hnd a model close to the observed variable. iLi & Oia^ 
(120 12h analyzed the photometry of V239, an Algol-type binary 
wit h P - 1.19 d, who se light curves were obtained earlier 
by IKaluznv et al.l (l2004l) . The solution indicates that the orbit 
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Table 7. BS mergers formed between 10 and 13 Gyr, evolved to the age of 13 Gyr 


model 

fmerg(years) 

M(Mq) 

log(L/Lo) 

logT, 

model 

fmcrg (years) 

M(Mq) 

log(L/Lo) 

logU 

0.75+0.4(1.5) 

1.037E+10 

0.981 

0.533 

3.849 

0.8+0.4(1.9) 

1.225E+10 

1.010 

0.479 

3.855 

0.8+0.4(1.7) 

1.040E+10 

1.024 

0.793 

3.842 

0.8+0.5(1.8) 

1.286E+10 

1.096 

0.635 

3.883 

0.85+0.5(1.8) 

1.087E+10 

1.153 

0.899 

3.785 

0.8+0.6(1.7) 

1.260E+10 

1.185 

0.872 

3.905 

0.75+0.4(1.6) 

1.146E+10 

0.974 

0.415 

3.844 

0.8+0.7(1.7) 

1.229E+10 

1.272 

1.506 

3.726 

0.75+0.5(1.5) 

1.143E+10 

1.065 

0.676 

3.875 

0.81+0.4(2.0) 

1.274E+10 

1.016 

0.460 

3.855 

0.75+0.6(1.5) 

1.194E+10 

1.148 

0.896 

3.893 

0.81+0.5(1.8) 

1.244E+10 

1.107 

0.731 

3.887 

0.8+0.4(1.8) 

1.133E+10 

1.016 

0.600 

3.860 

0.81+0.6(1.7) 

1.224E+10 

1.196 

1.082 

3.844 

0.8+0.5(1.7) 

1.178E+10 

1.105 

0.847 

3.879 

0.82+0.4(2.0) 

1.231E+10 

1.026 

0.532 

3.861 

0.81+0.4(1.9) 

1.183E+10 

1.021 

0.560 

3.860 

0.82+0.5(1.8) 

1.203E+10 

1.119 

0.867 

3.883 

0.81+0.5(1.7) 

1.143E+10 

1.116 

0.906 

3.778 

0.82+0.6(1.8) 

1.292E+10 

1.200 

0.834 

3.905 

0.82+0.4(1.9) 

1.141E+10 

1.034 

0.662 

3.864 

0.82+0.7(1.8) 

1.278E+10 

1.285 

1.016 

3.910 

0.83+0.4(1.9) 

1.103E+10 

1.045 

0.795 

3.856 

0.83+0.4(2.1) 

1.271E+10 

1.032 

0.517 

3.862 

0.83+0.4(2.0) 

1.187E+10 

1.038 

0.624 

3.866 

0.83+0.5(1.9) 

1.260E+10 

1.123 

0.754 

3.890 

0.83+0.5(1.8) 

1.167E+10 

1.130 

1.304 

3.749 

0.83+0.6(1.8) 

1.255E+10 

1.209 

1.013 

3.904 

0.84+0.4(2.0) 

1.144E+10 

1.050 

0.746 

3.866 

0.83+0.7(1.8) 

1.247E+10 

1.297 

1.309 

3.731 

0.85+0.4(2.1) 

1.180E+10 

1.056 

0.714 

3.871 

0.84+0.4(2.1) 

1.223E+10 

1.044 

0.599 

3.868 

0.85+0.5(1.9) 

1.177E+10 

1.146 

0.923 

3.804 

0.84+0.5(1.9) 

1.219E+10 

1.134 

0.903 

3.884 

0.86+0.4(2.1) 

1.127E+10 

1.068 

0.902 

3.831 

0.85+0.4(2.2) 

1.260E+10 

1.049 

0.585 

3.869 

0.86+0.5(1.9) 

1.130E+10 

1.158 

0.885 

3.831 

0.85+0.5(2.0) 

1.265E+10 

1.137 

0.798 

3.891 

0.865+0.4(2.2) 

1.178E+10 

1.069 

0.783 

3.871 

0.85+0.6(1.9) 

1.274E+10 

1.225 

0.981 

3.911 

0.865+0.5(1.9) 

1.105E+10 

1.164 

0.775 

3.870 

0.85+0.7(1.9) 

1.281E+10 

1.309 

1.062 

3.925 

0.865+0.5(2.0) 

1.189E+10 

1.157 

0.949 

3.817 

0.86+0.4(2.2) 

1.208E+10 

1.062 

0.698 

3.874 

0.8675+0.4(2.2) 

1.158E+10 

1.070 

1.835 

3.704 

0.86+0.4(2.3) 

1.282E+10 

1.056 

0.581 

3.870 

0.8675+0.5(2.0) 

1.158E+10 

1.161 

1.728 

3.708 

0.86+0.5(2.0) 

1.215E+10 

1.152 

1.016 

3.828 

0.87+0.4(2.1) 

1.148E+10 

1.068 

1.811 

3.705 

0.86+0.5(2.1) 

1.300E+10 

1.143 

0.758 

3.896 

0.87+0.4(2.2) 

1.179E+10 

1.055 

1.656 

3.710 

0.86+0.6(1.9) 

1.230E+10 

1.236 

1.540 

3.720 

0.88+0.4(2.2) 

1.151E+10 

1.069 

1.823 

3.704 

0.86+0.7(1.9) 

1.244E+10 

1.321 

2.126 

3.700 

0.89+0.4(2.2) 

1.123E+10 

1.088 

2.048 

3.696 

0.865+0.4(2.3) 

1.254E+10 

1.062 

0.638 

3.874 

0.7+0.4(1.5) 

1.226E+10 

0.927 

0.246 

3.829 

0.865+0.5(2.1) 

1.271E+10 

1.151 

0.834 

3.895 

0.75+0.4(1.7) 

1.261E+10 

0.968 

0.316 

3.839 

0.865+0.6(2.0) 

1.297E+10 

1.233 

0.965 

3.918 

0.75+0.5(1.6) 

1.269E+10 

1.056 

0.504 

3.866 

0.8675+0.4(2.3) 

1.276E+10 

1.057 

0.916 

3.807 

0.75+0.7(1.6) 

1.252E+10 

1.230 

0.919 

3.912 

0.8675+0.6(1.9) 

1.254E+10 

1.233 

0.881 

3.806 


Table 8 . Comparison of the observed parameters of binary V60 from M55 with one of CCBMs 


Parameter 

P(d) 

a{RQ] 

Ml (Mo) 

M2(Mo) 

MRq) 

Ri{Rq) 

log(L/Lo) 

logU 

V60 

1.18 

5.49+.05 

0.33+.02 

1.26+.03 

1.48+.01 

1.10+.02 

0.813+.013 

3.832+.008 

model 0.89+0.8(1.9) 

1.18 

5.41 

0.28 

1.26 

1.42 

1.13 

0.880 

3.854 


inclination is close to 90°, hence the resulting parameters are 
probably almost correct. The resulting component masses are 
M\ - 0.07 and M 2 - 1.20Mq. The model 0.8675+0.6(2.0) 
evolves after 13.5 Gyr into a binary with a period of 1.36 d 
and component masses 0.10+1.20 Mq. We stopped the calcu¬ 
lations at this age, but the mass transfer was still ongoing and 
the period was decreasing, so after the next few time steps the 
binary would resemble V239 even more. V228 in 47 Tuc is 
another Algol-type star. Its present parameters are P - 1 .15 d. 
Ml - 0.20 and M? = 1.51Mq (iKaluznv et al.l[200^ . ISarn^ 
(l2008h modeled this star using his evolutionary program. The 
program did not include magnetized winds but, instead, it con¬ 
tained two free parameters individually adjusted and describing 
the total mass and AM lost during the evolution. He concluded 
that the progenitor binary had a period of 1.35 d and compo¬ 
nent masses equal to 0.88 and 0.85 Mq. Allowing for mass loss 
by the winds and the present total mass of V228, we estimate 
that the initial component masses must have been close to 0.9 
Mq. None of the CCBMs is that massive. The metallicity of 47 
Tuc is also several times higher than used in CCBMs. However, 
the model 0.96+0.93(2.3) with the solar composition, calculated 
for another purpose some time ago, evolves after 12.4 Gyr into 
the configuration with period of 1.13 d and component masses 


0.20+1.49 Mq. This is quite close to the observed parameters of 
V228. 

V47 and V53 are CBs with short periods of 0.27 and 0.31 
d. The photometric solutions resulted in component masses of 
0.20+ 1.66 Mq for V47 and 0.74+0.91 Mq for V53 (iLiu et al.l 
I 2 OI ih . Unfortunately, the orbit inclination resulting from the 
photometric solution is equal to about 72°and 40°for V47 and 
V53, respectively, so quite far from 90°. In such cases the mass 
ratio is poorly constrai ned and can be far from the spectroscopic 
value (lRucinskill2006l) . The binary V47, with a total mass of 
1.86 Mq, is substantially more massive than any of the CCBMs 
discussed in this paper. In fact, the mass of t he primary is 
higher t han twice the TO mass (equal to 0.8 Mq. IKaluznv et aH 
(l2013ali '). which casts some doubt on its determination, although 
another route of producing a massiv e merger is possible, in¬ 
volving multiple mergers or collisions (Ivan den Berg et al.l200ll; 
ISandauist et al.ll2003t iLeigh & Sills|l201 lli . Clearly, a high mass 

of the primary is also difficult to reconcile with its low lumi¬ 
nosity Lpr = 2Lq (iLiu et al.ll201 ll) . The expected initial mass of 
another binary, V53, is about 1.8 Mq, which is again beyond the 
limits considered by us. Still, we can find models in the CCBM 
set with the period and q similar to V53, for instance, the model 
0.8675+0.8(1.9), which evolves after 11 Gyr into a binary with 


14 
































K. Stfpien and M. Kiraga: Blue stragglers and W UMa-type stars in globular clusters 


a period of 0.305 d and component masses of 0.66+0.86 Mq. 
We did not attemp t to find possible p r ogenitors of the bina¬ 
ries investigated by Kallrath et alJ d 19921) . [Helt et al.l ([T993h . and 
iRubenstein & Bailvnl (119961) because their parameters are even 
more uncertain. 

A more accurate search for progenitors of individual binaries 
in GCs will be the subject of a separate paper. 


4.2. Blue stragglers 

Recent progress in understanding the origin and properties of 
BS in GCs has been enormous, mostly du e to observa t ions o f 
GC cores w i th the Hubble telesc o pe of iPiott o et al. (f2004h. 


(l2Q13ah . Dalessandro et al. (l2013bl) . and lSanna et al.l (l2014l) . A 


iLeighet^ ( 2007), BeccarietaLJ (1201 2|). Dalessan dro et al 


comprehensive review of ob s ervati ons and resulting conclusions 
can be found in iBoffin et al.l (l2014l) . 

To compare models with observations, we first need to spec¬ 
ify which members of a GC are classified as BS. There is no 
single generally accepted method. To avoid contamination from 
photometric errors of MS stars at the TO point and blends 
(particularly in a cluster core), selection criteria require that 
a BS candidate be brighter than TO by a predefined quan¬ 
tity. In addition, a BS should lie at the extension of the clus¬ 
ter MS, but the width of the BS band is not strictly deter¬ 
mined and varies among different authors. Some authors plot 
boxes in the color-magnitude diagram ( ^ CMD), s ee, for instance , 
iLeigh et ~^ (1201 lah.jDale ssandro et al.l(l2013a) . lLi et an(l2Q13h . 
and Fiorentino et al l (|2014|) , others select BS by analyzing in¬ 
dividual candidates and a pplying specific criteria (iMoretti et al.l 
l2008tlFerraro et al.l2009h . The resulting set of BS approximately 
extends over 1-2 mag in V and 0.4-0.5 in B -V. Most of the ob¬ 
served BSs lie along the extension of the MS, but some stars, 
lying already beyond the TAMS, where evolved single stars and 
binaries can be found, are often included in the BS samples. The 
BS region adopted in this paper is bound by the ZAMS, the cor¬ 
responding isochrone (11, 12, or 13 Gyr), and the ZAMS line 
shi fted by AL - +0 8 (F igs.|6]El). 

iFerraro et al.l ( l2009l) detected a gap dividing BS in the clus¬ 
ter M30 into two sequences: blue and red. They argued that 
the separation results from two different formation mechanisms: 
the blue sequence, lying close to the ZAMS, originated from 
collision-formed mergers, whereas the red sequence that con¬ 
tains BS formed through mass transfer in binaries. Similar se¬ 
quences, albeit with a less pronounced gap betw een them, were 
detected in NGC 362 dPalessandro et al.1l2013bh . The blue se¬ 
quence agrees with the expected position of mergers formed by 
a collision of two stars with masses close to half a TO mass. 
Such stars are hardly evolved before merging, hence a freshly 
formed merger is expected to lie close t o ZAMS. The lin e sepa¬ 
rating both sequences was adopted after lTian et al.l (l2006h . It lies 
approximately at the extension of the MS binary sequence, that 
is, the ZAMS line shifted upward by 0.75 mag, and reproduces 
the expected position of unevolved, equal-mass binaries. It fol¬ 
lows from this interpretation that the blue sequence consists of 
single unevolved stars, whereas the red sequence should consist 
of binaries. In addition, the distribution of BS with a clear gap 
between both sequences requires t hat they were all for med in 
a single recent burst, identified by IFerraro et al.l (l2009l) with a 
cluster core collapse. The subsequent evolution of blue BSs will 
smear the distribution over CMD, so the gap should disappear. In 
addition, if BSs are f ormed continuously in a cluster, none of the 
features observed by IFerrar o et alJ ( ^09l) and iDalessandro et al.l 
(l2013bl) should occur ( Li et al. 20131) . To see how the models 



logTgff 


Fig. 9. Collective HR diagram of all mergers (filled circles) and 
all binary BS (asterisks) from Figs.|38] together with the corre¬ 
sponding isochrones. Rapidly rotating mergers are additionally 
marked with squares. The position of the eclipsing binary V60 
from M55 and of the model 0.89+0.8(1.9) are plotted as a di¬ 
amond and a plus sign. The ZAM S and the blue b oundary of 
the binary BS region determined bv lTian et alJ (l2006l) are shown 
as a dotted and broken line. Long dashes mark the adopted red 
boundary of the BS region. 


are situated relative to the blue and red sequences seen in M30 
and NGC 362, we plot all mergers from Tables [215] and |7] as 
filled circles in Fig. 0 and all binary BSs as asterisks. A few 
binaries with 1 < F < 3 d are also added. We do not list nu¬ 
merical data for them because the accuracy of the CCBM calcu¬ 
lations decreases with increasing period of an evolutionary ad¬ 
vanced binary, therefore their parameters are rather uncertain. 
The isochrones at 11, 12 and 13 Gyr are also plotted as solid 
lines. In addition, the ZAMS lin e (dotted), t he upp er edge of the 
binary BS region, determined bv lTian et al.l (l2006l) (broken line), 
and the adopted red boundary of the BS region (long dashes) are 
shown. The broken line delimits the region of binary BSs from 
the blue side. All CCBMs but one lie to the right of this line. The 
exception is model 0.8675+0.6(1.9). It became a CB at the age 
of 11 Gyr, but was lying below TO at that time (Table |2). By a 
continuing mass transfer, the gainer reached a mass of 1.161 Mq 
at 12 Gyr and the binary moved to the blue part of the BS region 
(Table |4]i. The gainer is hardly evolved at that age and the influ¬ 
ence of the redder secondary is insufficient to move the binary to 
the red part of the BS region. This indicates that in exceptional 
cases CBs can also be fo und in the blue sequence, as is observed 
in M30 and NGC 362 (IFerraro et al.l l2009t iDalessandro et al.l 
1201 3bl) . 

In contrast to binaries, the model mergers can be found 
in both parts of the BS region (Fig|9]l. Those with the lowest 
masses, around 1 Mq, lie close to ZAMS, along a narrow blue 
sequence. They were formed by two unevolved stars, hence re¬ 
sulting mergers are also unevolved even at the age of 13 Gyr. 
Their masses are between 0.9-1.1 Mq. A gap separating them 
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Fig. 10. Mass distribution of the BS mergers from Table|2](solid 
line) compared to the observed distribution of the pulsation 
masses of SX Phe stars in NGC 6541 (broken line). Both dis¬ 
tributions are normalized to unity. 


from the binaries is clearly visible. However, mergers more mas¬ 
sive than about 1.1 Mq are formed by two stars of which at least 
one has substantially depleted core hydrogen, which means that 
they are evolutionary advanced already immediately after for¬ 
mation. In effect, they arrive to the sequence tilted to ZAMS 
and extending to the red part of the BS region. Similarly, merg¬ 
ers formed collisionally in the evoluti onary models of a GC 
are scattered over the whole BS region (ISills et alj|2009l [20131 
IChatteriee et al.ll2013h . We conclude that binaries that move to 
the BS region as a result of mass ex change, arrive pre dominantly 
to the red of the line determined bv lTian et al.l ( l2006l) . However, 
the existence of a gap delimiting the narrow sequence of BSs 
close to ZAMS and extending to the most massive BSs cannot 
be reconciled with a continuous merger production. A burst in 
their product ion seems to be the most plausible explanation, as 
suggested bv iFerraro et al.l (l2009h . although this picture per se is 
not proof of the collisional formation of blue BSs. 

Masses of a few binary BS are known from the solution of 
the eclipsing light curves, sometimes supplemented with radial 
velocity curves (see previous subsection). Less is known about 
masses of single BS. Until recently, two methods have been ap¬ 
plied to determine their masses. The first was based on spec¬ 
tral and photometric obser vations, the second on the evolution¬ 
ary tr acks in the CMD ( Shara et al. 1997t iGeller & Mathied 
201 IL I2Q12I: iDe Marco et al.l 120051: iKaluznv et al.l 120071 12009 


_ Unfortunately, discrepancies o ften occur when both 

methods are used on th e same object (iDe Marco et"^ l2005t 
iGeller & Mathieu| j20j_2[). Recen tly, a new method was intro- 
duced by iFiorentino et al.l d2014h . It is based on the comparison 
of the pulsation properties of BS that belong to SX Phe - type 
variables. Applying this method to SX Phe stars in NGC 6541 
resulted in the pulsation masses of nine variable BSs. Six of them 
are pulsating in the fundamental mode and three in the first over¬ 
tone. For the last three stars the authors also calculated masses 


assuming fundamental pulsation, but they clearly preferred the 
first overtone identification. The masses agree very well with 
evolutionary masses, particularly when they are binned accord¬ 
ing to luminosity. We compare them with masses of the merg¬ 
ers from the BS region in Fig. |8] There are 45 such mergers 
in Table Q Figure [10] shows the mass distribution of the model 
mergers (solid line) together with the observed distribution (bro¬ 
ken line). Both distributions are normalized to unity. The agree¬ 
ment is surprisingly good, taking into account the low number of 
the observed masses and the fact that SX Phe instability strip is 
narrower than the BS region. The lack of the observed masses in 
the last bin, between 1.2-1.3 Mq , most likely results from this 
restriction. 


Recently, spectroscopic studies of BSs in a number of GCs 
have been performe d with the aim to investigate the m echanisms 
of their formation (iLovisi et al.l 1201iLovisil 1201 4 and refer¬ 
ences therein). One of the goals was to measure rotation rates 
of BSs. Freshly formed binary mergers are expected to rotate 
rapidly even though a g reat maiority of the orb ital AM is ex¬ 
pelled with some matter dRasio & Shapiroll995h . Rapid rotation 
is also expected for c ollision-formed mergers unless the colli¬ 
sion is strictly central (ISills et al.ll997il2()()ll) . But cool stars lose 
AM by the magnetized wind. In effect, they spin down from an 
initial rotation period of a fraction of a day to about one week 
in 10* years, that is, in terms of the eq uatorial veloc ity, from 
Veq ~ 100-150 km/s to less than 10 km/s (lBarnesl2009h . Except, 
possibly, for the most massive mergers, all other BSs in GCs are 
magnetically active and spin down on a timescale of 10* years. 
In Fig. 0 we highlighted all mergers formed in the last 10* years 
before 11, 12, and 13 Gyr with squares. Rotation evolution is not 
included in the CCBM code, but we assume that freshly formed 
mergers rotate rapidly. As we see, all these mergers lie close to 
the high-temperature boundary of the BS region, indicating that 
they are still unevolved. There are only nine rapid rotators in 
the figure, that is, about 10% of all BS mergers. This is an upper 
limit for the expected fraction of rapid rotators in a GC, achieved 
only if mergers formed in the last Gyr exist in the BS region. All 
mergers formed earlier will add to slowly rotating, thus decreas¬ 
ing the fraction of rapid rotators. For example, only three such 
rotators are expected to occur among 48 BS mergers in Fig. |8] 
(see Table |7]i, which is equivalent to about 6%. 


So far, accurate observations of vsin/ of several BSs in 


Lovisi et al. 

20 iqh 

Lovisi et al. 

201 3h 


|), NGC 6397 (iLovisi et al.l 120121) , and ’m 30 
The total number of BSs with a measured 
rotation rate in these clusters is 96, ten of which are rapid rota¬ 
tors with V sin i > 50 km/s (we exclude one W UMa-type binary 
counted as a rapid rotator in M30), this is about 10%. This is 
close to the predicted value (see above). Moreover, all rapid ro¬ 
tators in M4 and moderate rotators (with vsini > 10 km/s) in 
47 Tuc lie near the blue bou ndary of the BS region, indicating 
that th ey are freshly formed dFerraro et al.ll200ra IKaluznv et al.l 
l2013al) . However, this statistics is poorer when individual clus¬ 
ters are considered: eight rapid rotators are present in M4 (40% 
of the total measured sample in this cluster), whereas only two 
are members of three other clusters, of more than 80 stars mea¬ 
sured. This shows again that GCs are very different from each 
other, and a comparison of an individual cluster with theoretical 
predictions may be misleading. 


If mergers spin down as a result of magnetized winds, we 
expect a gap in the rotation distribution between rapidly and 
slowly rotating stars, resulti ng from the existence of the so- 
called Vaughan-Preston gap (lBarnesll201 ll) . The rotation distri- 
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butions of BSs show a gap (see Fig. 10 in iLovisi et al.ll2()T^ . 
although its statistical significance is low. 

Another goal of the spectroscopic survey of BSs was to 
search for anomalous surface abundances in BSs t hat would in¬ 
dicate the material proce s sed by the C NO cycle (iFerraro et alJ 
l2006t iLovisi et'aDl2010i 1201 2[ l2013h . As a star evolves, the 
abundance of carbon and decreases, whereas the abun¬ 
dance of nitrog en incre ases in the layers where nuclear re¬ 
actions operate (Ilbenlll9^ . The transition between the cosmic 
and modified abundance moves with time slowly outwar d un¬ 
til it r eaches the mass fraction ~ 0.3 for solar-type stars dlbenl 
Il967h . If enough matter is transferred to the companion, the ma¬ 
terial processed by nuclear reactions is exposed, hence anoma¬ 
lous abundances of these elements are expected to be present on 
the surfaces of both stars. Detection of any of these anomalies in 
a BS can be treated as an evidence for the mass transfer channel 
of its formation. An und erabundance of carbon ha s indeed been 
observed in field Algols (ISarna & De Grev3ll996ll . 

It follows from these considerations that at least 70% of the 
loser mass needs t o be transferred to its companion to e xpose al¬ 
tered abundances (llbenlll967l : ISarna & De Grevelll996h . In most 
CCBMs, only about 40-50% of the mass is transferred to the 
gainer before the two components merge. We should therefore 
not expect to see any CNO peculiarities in these binaries and 
mergers. A sufficient mass fraction is transferred only in a few 
binaries that ultimately merge as a result of the Darwin insta¬ 
bility. There are about ten such BSs among the binaries plotted 
in Fig. 0 All except one lie in the red part of the BS region 
or in the subgiant region. To these, one merger resulting from 
the model 0.8675-1-0.4(2.3) can be added because it was formed 
by the Darwin instability just before 13 Gyr. Its surface layers 
may also show carbon and oxygen depletion. Altogether, less 
than 10% of the plotted BSs are expected to show an underabun¬ 
dance of carbon and oxygen, which means that the probability 
of finding such an object among a random sample of BSs does 
not exceed 0.1. 

The observation s of GCs agree with th ese e xpectat i ons. T hey 
are summarized by IFerraro et al.l (1201 4l) and iLovis 1 IMIl- c 
and O depletion for a few BSs in 47 Tuc, M30 and w Cen was 
detected, but none of the observed BSs in three other clusters 
showed the expected alterations. The authors estimated that only 
about 10% of BSs may show CO depletion. That was fewer than 
they expected, therefore they speculated that the low percentage 
may result from the transient character of the alterations. While 
this assumption may still be viable, it is not needed to explain 
the observations. 


5. Discussion and conciusions 


We computed an extended set of CCBMs with initial param¬ 
eters expected to reproduce observations of CBs, NCBs, and 
BSs in GCs. The initial component masses were between 0.7- 
0.9 Mq for the primary and 0.4-0.8Mq for the secondary 
components. The initial orbital period distribution had a short- 
period cut-off at 1 .5 d resulting from the accepted binary forma¬ 
tion mechanism dBos ^ 119931: lBonnelllll994 iKratter et al.ll20Tol: 
iMachida et al.l 12008 ) and a maximum between 2-3 d resulting 
from the KCTF mechanism operating on suitably shaped triple s 

dEggleton & Kiseleva-Eggletor^ 200^|Perets & Fabrvckvl2009ll . 

The metal abundance was fixed at Z=0.001. Magnetized winds 
from both components were adopted as the driving force for or¬ 
bit evolution. Evolutionary calculations were carried out to the 
age of 13.5 Gyr, but detailed analysis of the evolved objects was 
performed for three ages: 11, 12, and 13 Gyr. 


The results show that all binaries with primary masses lower 
than 0.7Mq and almost all with masses 0.75 Mq remain de¬ 
tached untill the age of 13.5 Gyr. On the other hand, almost all 
binaries with initial primary masses close or equal to 0.9 Mq 
complete the evolution as mergers or common-envelope prod¬ 
ucts before they reach the age of 11 Gyr. Only binaries with 
primary masses embracing the TO mass at a specified age of 
11, 12 or 13 Gyr and with initial periods from a narrow inter¬ 
val around 2 d form CBs or NCBs with periods shorter than 1 
d, which should be observed as W UMa-type stars in a GC of 
that age. Most of these variables do not survive for more than 
1 Gyr and either form a merger or an Algol-type star with a 
period of several days. Consequently, only 4-5% of the mod¬ 
eled binaries are in contact or near-contact at each specified age. 
Observations show, i ndeed, that W U Ma-type binaries are rare 
among GC members (lRucinskill2000l) . Moreover, about 90% of 
model CBs lie close to, but below TO point on the HRD, which 
means that on average, only one in ten CBs falls in the BS re¬ 
gion. Considering the faint limit for CBs, the absence of model 
CBs fainter than ~ 1 mag below TO point agree s well with ob¬ 
serva t ions of GCs (iKaluznv & Thompsonll200^ iKaluznv et al.l 
1201 0112013ah . 


As the calculations show, freshly formed mergers have 
masses between 0.9-1.5 Mq, as compared to a typical TO 
mass of ~ 0.8Mq in a GC. Thus, all they arrive in the BS 
region. In effect, they outnumber binary BSs by at least an 
order of magnitude, so among the random sample of BSs, 
fewer than 10% are expected to be close binaries. The observa¬ 
tions support this result, although substanti al fluctuations of this 
fraction occur among individual cluster s ( Eerraro et alJ 1200^ 
iDalessandro et alJ 1201 3bl IKaluznv et alJ l2010l l2013ah . Algols 
with periods longer than 1 d occur predominantly to the red of 
the BS region because of the influence of a cool, Roche-lobe 
filling component that reddens the binary. 

All mergers with masses lower than 1.1 Mq are still un¬ 
evolved and lie along ZAMS, thus forming a blue BS sequence. 
On the other hand, all binary BSs lie in the red part of the BS 
region, thus forming a red sequence. A gap separates the two 
sequences. The existence of such a gap has been noted in two 
GCs, but it was linked to a highly nonuniform B S formation rate 
(lEerraro et al.ll2009t [Dalessandro et alj|2013bh . Eurther investi¬ 
gations should resolve this discrepancy. 

The uncertainties of the CCBMs mostly stem from the un¬ 
certainty of the coefficients in Eqs. |2ll8]and the accuracy of the 
stellar models. An increase of the AML rate by 30% results in a 
decrease of the orbital period at RLOE by about 0.2-0.3 d. This 
has a different influence on binaries with RLOE occurring on 
the MS and occurring beyond the TAMS. In the former, RLOE 
takes place earlier, when the primary is less evolved, hence the 
mass ratio past rapid mass exchange is higher by about 0.05-0.1. 
Component merging and its following evolution also occur ear¬ 
lier. In effect, CBs observed at the specified age originate from 
binaries with initial periods longer by the same amount. In the 
latter, the shorter period at RLOE has a negligible effect on the 
process of mass exchange because the fast radius increase of the 
subgiant component changes the time of RLOE by a mere few 
time steps. The mass ratio past rapid mass exchange also changes 
little, so only the orbital period is shorter, again, by the same 
amount. A decrease of the AML rate by 30% has the opposite 
effect to the same extent on the binary parameters. In short, the 
uncertainty in the AML rate of 30% translates into an uncertainty 
of about 10% in the binary and merger parameters. 

A decrease of mass-loss rate by a factor of two has an even 
weaker effect. The total mass lost of a CB is between 6-7% of 
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the initial binary mass. Its cut in half shortens stellar evolution¬ 
ary time scales by about 10% so RLOF occurs somewhat earlier, 
but the hnal parameters of the considered models agree within 4- 
5% with those presented above. An increase of the mass-loss rate 
by a factor of two influences the resulting parameters stronger, 
but, again, the change is lower than 10%. Only a substantially 
higher mass-loss rate (e. g., by a factor of 5) would qualitatively 
change the developed picture. Spherically symmetric mass loss 
always causes a period increase. For each AML rate there exists 
a maximum value of the mass-loss rate such that for higher mass- 
loss rates orbit widens. That would substantially change the bi¬ 
nary evolution. There are no accurate measurements of mass- 
loss rates for low-mass stars, therefore an extrapolation of data 
for solar-type stars was used. The mass-loss rates from Eq. [T] 
give between 0.5- 1.5x10“"Mq per year for a 0.8 Mq star at the 
ZAMS and at the TAMS, respectively, and 2-7xlO“'^M0 per 
year for a 0.5 Mq star. Other a uthors suggested s imilar rates for 
low-mass dwarfs, for example. ICranmer &_Saai ( 201 Ih. or even 
some what lower (~ IO^^^Mq per year. Ivan den Oord & Dovlel 
Il997h . It seems therefore that a substantially higher mass-loss 
rate is unlikely. 

An assessment of the uncertainties connected with stellar 
evolutionary models is more difficult. Even the models of sin¬ 
gle stars differ significan tly among different authors because 
of diff erent input physics dVerbunt & Phinne~^ll99-5tlNataf et alJ 
l2012h . In particular, the timescales of the consecutive evolution¬ 
ary phases and radii of the models may differ by several percent. 
Of the same order is the uncertainty in stellar radii, resulting 
from magnetic infl ation observed in rapidly rotating cool dwarfs 
(iTorres et alj|2010l) . Eortunately, this uncertainty has only a mi¬ 
nor effect on the final results. The age of RLOE and the merging 
time change by several time steps when we change the stellar ra¬ 
dius by 10%, and the period boundary between CBs and NCBs 
(which,anyway, we treat together) shifts only little. 

These considerations about binary BSs only relate to close 
binaries. Any BS can also be a member of a long-period binary. 
In fact, our assumptions require that most of these objects pos¬ 
sess (or possessed in the early evolutionary phases) a tertiary 
component that shortens the inner binary period by the KCTE 
mechanism, as suggested by IPerets & Eabrvck\l (l2009l) . These 
authors also discussed the observational consequences of this as¬ 
sumption. They compared the timescales for dynamical disrup¬ 
tion of triples in cores and outer regions of stellar clusters with 
the timescale for a formation of a close binary. Both time scales 
are functions of the outer binary period. It follows from the com¬ 
parison that for Pout < lO'* d dynamical encounters are ineffec¬ 
tive in the outer parts of a GC. Thus a large part of BSs should 
occur there in triples (before merging) or in binaries. Only those 
formed in initially soft triples might be disrupted. Altogether, the 
fraction of binaries and triples among BS should be higher than 
the overall binary fraction in the environment where they are ob¬ 
served, and most should be hard. The much higher frequency 
of the dynamical interactions in GC cores make the prediction 
more difficult. The timescale for the close binary production is 
short enough for Pout < 10"^ d, so the KCTE mechanism can 
still operate in a GC core, but a high interaction rate results in 
a disruption, a change of the conhguration or a merger of the 
inner binary. The outer binaries are even more susceptible to 
disruption or component exchange. In effect, the fraction of bi¬ 
nary BSs in a GC core is expected to be lower and they will 
have different parameter distribution from those in low-density 
regions. If the triple system is not disrupted, a second-generation 
BS may be formed when a hrst-generation BS, considered here, 
transfers matter to the original tertiary after reaching the red. 


or asymptotic giant branch. In effect, a population of BSs with 
white dwarf companions is expected. The observations of the 
old open cluster NGC 188 revealed that a signih cant fraction of 
BSs has such companions (IGosnell et al.ll2014li . Unfortunately, 
no such observations for GCs are available. 

Eor very few individual CBs in GCs, sufficiently accurate 
parameters are known. Although we did not attempt at present to 
hnd progenitors that fully match the observations, we can always 
hnd a model within the CCBM set that closely reproduces the 
observed parameters. Eor example, the model with initial masses 
0.89-1-0.8 Mq and the initial orbital period of 1.9 d evolves after 
11 Gyr into a binary with masses 1.26-1-0.28 Mq and a period of 
1.18 d, which matches the variable V60 from M4 very well, for 
which the most accurate data have been determined (see Table [8] 
and Eig.|9]l. Adding age and metallicity as additional constraints 
should result in a unique initial model for this binary. 

One of our main results shows that the binary-originated BS 
cover most of the region on the HRD where the observed BSs 
occur. The least massive mergers lie close to ZAMS. This is due 
to the low masses of the merging components with hardly de¬ 
pleted core hydrogen. The product of merging is then unevolved 
and can remain in the BS region for 2-3 Gyr before moving 
to the red giant region. They are indistinguishable in the HRD 
from the collision-originated mergers of two low-mass stars. The 
higher the component masses, the more advanced evolutionally 
the merger, hence the blue boundary of the BS sequence turns 
more and more away from ZAMS (see Eigs. IMSl) . Evolutionary 
advanced mergers with core hydrogen depleted and binaries past 
RLOE populate the red part of the BS region. Only about 6-10% 
of the mergers are expected to rotate rapidly, and all probably 
lie close to the blue BS boundary. Similarly, a low fraction of 
BSs may show chemical alterations resulting from exposing lay¬ 
ers where the CNO cycle operated. They should occur predomi¬ 
nantly among binaries with low mass ratios where a sufficiently 
large amount of mass has been transferred between the compo- 
nents. These predictions agree with observa t ions o f a few GCs 
(lEerraro et al .1120061 ll .ovi si et al .1120 1 (112011120 1 3h . 

To conclude, none of the observed parameters of an individ¬ 
ual single BSs can uniquely point to its origin: binary or col- 
lisional. The presence of binaries among BSs obviously indi¬ 
cates that at least a fraction of these objects has a binary ori¬ 
gin. Can this be true for most BSs? We noted profound differ¬ 
ences among individual GCs regarding the properties of BSs, 
therefore a single picture that would be applicable for all clus¬ 
ters seems improbable. Evolutionary modeling of a GC, includ¬ 
ing stellar and binary evolution, should provide credible predic¬ 
tions ahouUhejelStiveJmpqrtanceofdifferent_groduction chan¬ 
nels (iHvpki & Gierszl2013HChatteriee et al.l2013h . On the other 
hand, a stati stical analysis of GCs should help to verify these 
predictions (iMilone et al.l I2012I : ISills et aP l2013h . It seems at 
present that models of GCs favor for mation mechanisms involv¬ 
ing binaries. iHvpki & Giersd (l2013h modeled a cluster with 10^ 
objects, 20% of which are binaries. They considered ten different 
ways of BS formation. The results showed that single-star colli¬ 
sions produce a negligible number of BS, 13.6% are formed by 
collisions involving binaries (single-binary and binary-binary), 
and the rest by binary mergers. They noted, however, that these 
results may not apply to more massive and dens er clusters. 

So mewhat different results were obtained bv IChatteriee et all 
(12013h . who showed that in a dense environment (above 
IO^Mq/pc^) most of the BS are formed by collisions mediated 
by binaries, but in lower-density clusters most BS (up to 60%) 
result from mass transfer in binaries. In effect, we can expect 
a varying proportion between both kinds of BSs among differ- 
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ent clusters. This is in line with iFerraro et alJ (l2Q12h . who ar¬ 
gued that BSs formed in cluster cores originate from collisions, 
whereas those in cluster peripheries result from binary evolution. 
Therefore a flat density distribution is seen in dynamically young 
clusters, whereas an increasing concentration toward the cluster 
center is visible in dynamically older clusters when more mas¬ 
sive BSs from the outer parts gradually migrate into the core. In 
addition. lLi et alJ (1201 3h announced the observations of bimodal 
distribution of BSs in the Large Magellanic Cloud GC Hodge 11. 
They noted that BSs from the central distribution peak lie within 
the blue sequence in the HRD, whereas those from the outer 
peak are placed in the red part of the BS region. Assuming that 
the position in the HRD uniquely determines an origin of a BS, 
they concluded that core BSs are collisionally formed, whereas 
binary-formed BSs dominate in the outer parts of the cluster. The 
results of the present paper show, however, that this assumption 
is not correct. We also note that the statistical comparison of 
the BS population in several GCs with the simplihed formation 
model clearly indicates that bina ry-formed BSs dom inate those 
collisionally formed everywhere (iLeigh et al.ll201 Ibl) . 

In brief, the present evolutionary model of cool close bina¬ 
ries fully explains the origin, evolutionary stage, statistical prop¬ 
erties, and individual physical parameters of CBs and NCBs ob¬ 
served in GCs. 

Our results favor a binary origin of BSs, but this is still an 
open question. More observations of these stars in different re¬ 
gions of GCs, in particular determinations of masses and multi¬ 
plicity fraction compared to the environment, will shed light on 
their origin. Better evolutionary GC models will help to interpret 
observations and to understand substantial differences observed 
among different GCs. 
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